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ABSTRACT 

We explore brane induced gravity on a 3-brane in six locally flat dimensions. To regulate 
the short distance singularities in the brane core, we resolve the thin brane by a cylindrical 
4-brane, with the geometry of AD Minkowski x a circle, which has an axion flux to cancel 
the vacuum pressure in the compact direction. We discover a large diversity of possible 
solutions controlled by the axion flux, as governed by its boundary conditions. Hence brane 
induced gravity models really give rise to a landscape of vacua, at least semiclassically. For 
sub-critical tensions, the crossover scale, below which gravity may look AD, and the effective 
AD gravitational coupling are sensitive to vacuum energy. This shows how the vacuum 
energy problem manifests in brane induced gravity: instead of tuning the AD curvature, 
generically one must tune the crossover scale. On the other hand, in the near-critical limit, 
branes live inside very deep throats which efficiently compactify the angular dimension. In 
there, AD gravity first changes to 5D, and only later to QD. The crossover scale saturates at 
the gravitational see-saw scale, independent of the tension. Using the fields of static loops 
on a wrapped brane, we check the perturbative description of long range gravity below the 
crossover scale. In sub-critical cases the scalars are strongly coupled already at the crossover 
scale even in the vacuum, because the brane bending is turned on by the axion flux. Near 
the critical limit, linearized perturbation theory remains under control below the crossover 
scale, and we find that linearized gravity around the vacuum looks like a scalar-tensor theory. 
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1 Introduction 



1.1 Prologue 

The mystery of the cosmological constant is arguably the most pressing and puzzling problem 
in contemporary fundamental physics. To date, the attempts to explain a small cosmological 
constant in the framework of effective field theory formulation of matter coupled to gravity, 
based on conventional lore of naturalness alone, have not yielded an answer (see the classic 
work [1] which still provides a state-of-the-art review). On the other hand, cosmological 
observations [2] strongly suggesting a small cosmological constant, and the discovery of the 
landscape of string vacua [3], have lent support to the idea of statistical selection of vacua, 
and the value of the cosmological constant in them [4] -[9], prompting a new debate about 
anthropic reasoning in physics. 

One may try to alter the problem by embedding our universe in a fundamentally higher- 
dimensional space-time in some way [10]-[12]. If, for example, our universe is a brane in 
extra dimensions, it would also have extrinsic curvature, and then one may divert the brane 
vacuum energy into the extrinsic curvature [12]. In this way, the vacuum energy could remain 
invisible to long distance 4D gravity, yielding a very weakly curved universe. However, while 
attempting to recover 4D General Relativity at large distances by compactifying the bulk, 
with covariantly conserved sources, one needed additional branes in the bulk, that restore 
a fine-tuning similar to the standard 4D one [12, 13]. The resulting picture was again that 
of a landscape: a theory admitting multiple 4D vacua, classified by the free integration 
constants for some of the bulk fields (see, e.g.. [14]). Using these integration constants, one 
could hope to derive a framework where the effective cosmological constant could change 
in small amounts from one background to another. But one gains no more than that: the 
landscape reemerges, and one must retreat to either statistical or anthropic arguments to 
select a phenomenologically viable ground state. 

A different idea for recovering 4D General Relativity at large distances has been pursued 
in the so-called brane induced gravity theory [15]. In more interesting variants of brane- 
induced gravity, the extra-dimensional space has infinite volume. When the bulk volume is 
infinite, the 4D graviton zero mode is not normalizable, and hence it completely decouples. 
Thus 4D gravity ought to emerge from the exchange of the continuum of bulk modes, which 
are all massive from the 4D viewpoint. To ensure this, [15] introduced induced curvature 
terms on the brane, arguing that they will be generated by brane quantum corrections 
anyway. If the scale M 4 that normalizes them is much greater than the bulk Planck scale, 
the brane kinetic terms would pull the bulk gravitons with wavelengths shorter than a 
certain crossover scale r c very close to the brane, yielding the momentum transfer due to 
the scattering of virtual bulk gravitons oc 1/p 2 for the momenta p > r c _1 , and a very small 
coupling of the three-point vertex ~ 1/M 4 . This gives a force which scales as 1/r 2 , simulating 
AD behavior [15], and is nicely illustrated in the exact gravitational shock waves of [16]. 

This force contains admixtures of longitudinal, helicity-0 gravitons, and sometimes also 
radions. They couple gravitationally and their long range contributions to the force never 
completely disappears in perturbation theory [17]. However it has been argued, for massive 
gravity [18] and similarly for brane induced gravity [19], that the theory goes nonlinear 
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at large distances from the source, and that the nonlinearities might screen away the extra 
scalar modes. On the other hand, this also indicates [20]- [23] that brane induced gravity runs 
into a strong coupling regime at macroscopic distances r strong ~ (r^/A^) 1 / 3 , requiring very 
low scale UV completions, and forcing the question of whether such completions even exist! 
The extra scalars may also turn into a perturbative ghost on self-accelerating backgrounds 
in 5D bulks [22, 24, 25, 26, 27, 28]. This suggests that the self-accelerating backgrounds 
are unstable, but the strong coupling at subhorizon scales obstruct the exploration of their 
fate in perturbation theory. Still, the shock waves [16, 25] can identify energy leaks from a 
self-accelerating brane into the bulk, that render the dynamics manifestly different from a 
AD one. 

Similar mechanisms may also operate on higher-codimension defects [15, 29]. The authors 
of [29] have pursued ideas how UV regularizations of gravity might affect the crossover 
scale and the IR spectrum of the theory, hoping to suppress the extra graviton helicites. 
Subsequent discussions involved arguments as to how some regularizations may suffer from 
ghosts [30], which precisely cancel the extra graviton helicities, and how to avoid them 
[31, 32]. The explicit constructions are hindered by the short-distance singularities in the 
core of the defect, which must be regulated before one can reliably calculate the low energy 
behavior of gravity. 

We feel that exploring brane induced gravity on codimension-2 brane models is partic- 
ularly interesting. In this case we can find the exact background solutions that play the 
role of the AD vacua directly, thanks to the magic of gravity in 3D. It is well known from 
early braneworld analysis that 3-branes in 6D need not locally curve the bulk even if they 
carry tension [33]. Although this leads to flat AD geometries for nonvanishing tensions, it 
does not help with the cosmological constant problem after one compactifies the bulk to get 
long range AD gravity. This was immediately realized already in [33] to be in agreement 
with Weinberg's no-go theorem [1], and was discussed in more detail in [34]. Still, many 
interesting properties of gravity localized on codimension-2 defects, flat or not, in a 6D bulk, 
compact and not, have been since explored in [35]- [52], including exact black holes and gravi- 
tational shock waves straddling the brane [50]. The sheer proliferation of such configurations 
points to a structure as rich as the landscape of supergravity solutions, albeit more exotic. 
Sketching a chart of this realm is the primary aim of the present paper, as we now elaborate. 

1.2 The Landscape of Brane Induced Gravity 

We begin our foray into the landscape of brane induced gravity by exploring thin static 
codimension-2 branes, which exist for sub-critical tensions A < A c = 2nMg and have a flat 
AD induced metric. Because the brane induced curvature terms identically vanish, these 
solutions are identical to the ones in 6D General Relativity with an empty bulk. We then 
construct exact gravitational shock waves of relativistic particles on the brane, generalizing 
the solutions of [50] to the framework of brane induced gravity. The shocks confirm that 
brane induced gravity on singular defects remains AD all the way to infinity [15, 29]. 

Yet, these solutions are pathological. Even an infinitesimal displacement of the probe 
from the brane into the bulk reveals that the thin brane theory is singular, since the grav- 
itational field jumps from a finite value to exactly zero everywhere in the bulk [15, 29]. 
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Mathematically, the discontinuities on thin branes can be attributed to an ill-posed exterior 
boundary value problem: modelling brane sources by points on an infinite cone introduces 
singularities both near and far. One turns it into a well-posed problem simply by replac- 
ing the point sources with finite rings, and considering the combined interior and exterior 
problems. Physically, this amounts to smearing the source over a finite space. We do it 
by replacing the thin 3-brane by a 4-brane wrapped on a circle, by adding an axionic field 
on the brane, whose vev breaks translations along the circle [43, 51]. This is similar to the 
Scherk-Schwarz mechanism [53] for supersymmetry breaking and mass generation by dimen- 
sional reduction, an important ingredient of the string landscape model-building [3]. Much 
like in the string landscape, the axion 'charge' determines the radius of the compact circle 
r for a given value of the tension, and we must choose it carefully to get the desired value 
of ro- Alternatively, if some other dynamics were to fix ro, we need to tune the axion charge 
to precisely cancel the pressure of the brane tension along the circle, because otherwise the 
background could not simulate the AD Minkowski vacuum from the get go. 

The thick branes will still have a flat static geometry with an infinite bulk, but only for 
tensions smaller than the critical value. When the tension exceeds the critical value, our 
solutions show how the bulk compactifies and develops a naked singularity that soaks up the 
brane's gravitational field lines, just like the supercritical cosmic string in AD [54]. These 
solutions look like singular teardrop compactifications analyzed by Gell-Mann and Zwiebach 
in [55]. As with AD supercritical defects [54], one expects that for supercritical branes, as 
well as for the cases with a mismatch between the tension and the axion charge, there exist 
nonsingular solutions which describe some curved, nonstationary backgrounds. Examples 
which describe topologically inflating defects have been found recently in [56]. 

We next calculate the crossover scale where gravity changes from 4D, again using the 
shock wave solutions, sourced by relativistic strings moving on the wrapped brane. Below 
the crossover scale, the shock wave indeed approximates the AD Aichelburg-Sexl solution 
[57] down to r , whereas at much larger distances it changes over to the 6D shock of [58], for 
sub-critical branes. The physics of crossover is very intricate: it is controlled by the deficit 
angle of the vacuum solution. This is because of the 'lightning rod' amplification of the bulk 
gravitational coupling, leading to Mq cS = (1 — &)Mg [50], where 1 — 6=1 — measures 
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the deficit angle of the brane. It is intuitively clear why it affects the gravitational force: 
on a cone, gravitational field lines spread more slowly than on a plane, and hence gravity 
must look stronger [50]. For sub-critical strings with deficit angle smaller than 2-7T, our exact 
crossover scale is qualitatively very different from the see-saw scale of [29], yielding instead 
r 2 ~ J^ 4 , up to a log, where M 4 is the effective AD Planck scale found by compactifying 
the wrapped 4-brane theory on the circle of radius r . 

On the other hand, in the near-critical limit the bulk compactifies to a cylinder, which 
opens up into a cone only very far from the brane, at distances > j^. In the cylindrical 
throat gravity is in a 5D gravity regime, separating the AD and 6D ones. Thus the crossover 
scale beyond which gravity is not AD is really the demarcation between the AD regime and 
5D gravity which lives inside the throat [59]. To see where this happens, we can use the 
naive crossover formula on a codimension-1 brane, realized as a wrapped 4-brane in a 6D 
flat space with a compact circle: by Gauss formulas for Planck masses M\ e ^ = M|r and 
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M| e// = M 6 4 r , we find exactly r c ~ M t e// ~ jjb, which is our crossover formula from the 

Be// 6 

exact shock waves, and which is the same as the see-saw scale r c ~ Mi 4 of [29]. Thus we 
see precisely how the see-saw mechanism emerges. 

Our analysis reveals the presence of a large topographic diversity of solutions. The axion 
flux which cancels the pressure in the compact direction sets the radius of the circle by 
q 2 = 2A 5 rQ. Since the AD Planck scale is Mj ~ M|r ~ ^|=, and, for sub-critical branes, 

the crossover scale r c ~ ^j2\^M^-2\^q we can ^ n< ^ ^ e cnar S e Q that wm produce any desired 
M 4 for any values of M 5 and A 5 . But then, the crossover scale is uniquely fixed in terms of 
these parameters. The vacuum energy problem reappears in brane induced gravity: although 
the brane is flat, when we fix M 4 , we must finely tune both the axion charge and the tension 
of the brane to get r c to take some desired value on a sub-critical brane. On the other 
hand, once we fix A 5 , M 5 and M 6 , and let the charge q vary, we will find that the thickness, 
the effective 4D Planck scale, and, for sub-critical branes, the crossover scale of the brane 
also vary, spanning a wide range of theories below the crossover scale. The near-critical 
branes, where the deficit angle approaches 2it, are more interesting since their crossover 
scale saturates at r c ~ j^, which is completely independent of the brane tension in the 

leading order. This is because their conical throat [59] shields them from the asymptotic 
infinity. Only their effective 4D Planck mass is sensitive to the brane tension, and changes 
when the tension varies and the axion charge is fixed. The charge and the tension still need 
to be tuned right to get a desired value of M 4 , but this may be substantially easier. 

This is a signature of a landscape: low energy parameters of a theory depend on boundary 
conditions, instead of being uniquely fixed by a symmetry principle. We could also imagine a 
brane where q and A 5 might slowly vary along the brane, leading to regions of 4D space with 
a different strength of simulated 4D gravity, but always with a flat static 4D vacuum for 
sub-critical vacuum energies! Thus, we unveil a classical landscape of codimension-2 brane 
induced gravity, where the brane may remain flat since the curvature is soaked up by the 
deficit angle, but M 4 and generically r c depend on the vacuum energy. 

Including quantum mechanics in this picture requires a lot of care. An immediate, semi- 
classical issue is that since the axion is a phase of some complex 5D scalar field, its charge 
is quantized in the units of a new UV mass scale ji. The quantization law, q = /i 3 / 2 n, where 
n is an integer winding number, yields M 4 ~ ^= n, so that for a given set of dimensional 



parameters we must choose the right winding number. The value of M 4 changes discretely 
with it. Similarly, for sub-critical branes, the crossover scale can only change in discrete 
jumps on static backgrounds. Thus the static configurations will only exist at discrete loca- 
tions of the landscape, most of which is covered by nonstationary solutions, like in the string 
landscape with a discretuum of vacua [3]. Finding the nonstationary solutions is outside of 
the scope of this work; we hope to return to it elsewhere. Nevertheless as long as gravity 
is treated classically and field theory is consistently cut off at some scale in the UV, this 
can be done in principle, showing that the classical landscape which we glimpse at can be 
extended semiclassically. However, there remains the looming question about UV comple- 
tions of brane induced gravity models, and their embeddings into the theories of quantum 
gravity [15, 60, 61, 62]. We will have nothing new on this to add here, taking a more pedes- 
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trian approach of simply charting out the (semi) classical countenance of what the quantum 
landscape might be. 

The shock wave solutions have taught us that at distances between r and r c , the theory 
contains 4D General Relativity. To see what else is there, we develop, in painstaking detail, 
linearized perturbation theory about a wrapped brane. For sub-critical tension the scalars 
are strongly coupled essentially exactly at the crossover scale r c , which is the Vainshtein scale 
[18] of the vacuum itself. This happens because the scalar gravitons probe the asymmetric 
vacuum energy distribution, set up by the background axion flux, which triggers brane 
bending. However in the near-critical limit linearized perturbation theory around the vacuum 
remains under control below the crossover scale. Solving for the linearized gravitational field 
of a static uniform ring of mass on the brane, we find that below the crossover scale the 
theory contains helicity-0 and radion-like scalars in addition to the helicity-2 modes. The 
radion modes decouple but the helicity-0 modes remain active at the linearized level. Thus 
the linearized theory around near-critical vacua approximates a Brans-Dicke gravity with 
uj = 0. This would disagree with the classic tests of General Relativity [17], but perhaps 
non-linearities or dynamics on vacua of more complex structure could come to the rescue. 
We however find that rings of brane matter built of the lightest KK states on the wrapped 
brane do not entice any fast instabilities in the leading order of perturbation theory. While 
we suspect that fast instabilities might remain absent in perturbation theory beyond linear 
order, we have not proven it, and it remains a question for future studies. 

The paper is organized as follows: we begin with a lightning review of thin codimension- 
2 flat branes and their short distance singularities, as revealed by shock waves, in section 
2. We then construct the regulated backgrounds, with a 4-brane wrapped on a circle, in 
section 3. We also discuss how critical and supercritical branes compactify the bulk, the 
latter inducing a naked singularity far away. We derive regulated shock wave solutions in 
section 4, to explore the long range properties of gravity and compute the crossover scale. 
In section 5 we set up linear perturbation theory, and give the linearized solution of a static 
ring of mass on the 4-brane. We end with conclusions. 

2 Tensional Thin Branes on Bulk Cones 
2.1 Vacua 

The field equations describing a 3-brane in an empty 6D bulk in brane induced gravity are 
[15, 29], in a brane fixed Gaussian-normal gauge, 

M%G & A B + MlG^ u 5 A ,5 v B 5^\y) = T\5 A ^ B 5^\y) . (1) 

Our conventions are that the indices A, B, . . . count all spacetime dimensions, while fi, u, . . . 
run over the brane worldvolume. The 5-function is the tensor 8^ 2 \y) = : ^U5(yi), given by 
the normalization of the usual tensor density by the metric determinant of the extra space, 
coordinatized by (1/1,1/2) (we will use for the coordinates along the brane worldvolume). 
As in [15, 29] we take the stress-energy tensor T^ v to be localized on the brane, separating 
the tension A from the matter sources by T^ v = — A^ + r'V 
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The vacua of this theory are maximally symmetric solutions with r^ v = 0. Substituting 
this in (1), imposing a flat brane g^ v = r]^ and tracing over, we find R 2 = jjk5^(y), 

where R 2 is the curvature in the spatial dimensions transverse to the brane. The solution 
is a conical space [33], which can be seen most simply as follow. Since any metric in two 
dimensions is conformally flat, we can write ds 2 2 = e^^dy 2 . Then by conformal flatness of the 
metric, R 2 = e^V^i?, and recalling the Euclidean 2D Green's function \n(\y\/£), which obeys 

Vy\n(\y\/£) = 2irH5(yi), we see that to get the solution we can set -d = 2b\a.{\y\/i) to find 
R 2 = ATTbe 2bhl ^^m( yi ). We next compare this to R 2 = j^5^(y) = ^e 261n ^/^n5(^) 

6 6 

from above, from which we finally obtain 

(2) 



2vrM t 

The arbitrary length scale £, needed for dimensional reasons, completely drops out of these 
equations. Then substituting the solution for going to spherical polar coordinates and 
changing the radial coordinate to p = -^£ b \y\ l ~ h , we get 

ds 2 2 = dp 2 + (1 - b) 2 p 2 d<j) 2 . (3) 

Since the range of the angular variable is 2n, 1 — b is the deficit angle induced by the 
tension, as per Eq. (2). The full QD vacuum solution is 

dsl = rj^dx^dx" + dp 2 + (l- b) 2 p 2 d<p 2 . (4) 

This is identical to the thin 3-brane in a flat bulk found in models with large extra dimensions 
[33] because the induced curvature identically vanishes on flat branes. The bulk geometry 
is a cone, given in Fig. (1). It is well known, however, that for backgrounds with matter 
sources on a thin brane in higher- dimensional gravity the short distance structure of the 
geometry near the brane is more involved [46]- [48, 50, 51]. In general, the length scale I 
does not drop out but needs to be promoted into a UV regulator, smoothing the tip of the 
cone in Fig. (1). As we will see below this becomes even more important in brane induced 
gravity. 

Clearly, we see that the solution (4) is not well defined at b = 1, corresponding to the 
critical value of the tension A cr = 2ttMq. While the metric (4) makes sense, as it stands, 
for b > 1, it is not clear what it represents since we can't smoothly deform it to relate it to 
sub-critical cases. Hence in the critical limit, and beyond it, we cannot tell what happens 
to the brane until we sort out how to match the brane's core to the exterior, in precisely the 
same manner as for the critical and super-critical strings [54]. 



2.2 Shocks on Thin Branes 

A simple way to probe the nature of gravity, as argued in [16, 50], once a vacuum solution is 
known, is to shock it. By this we mean, put a relativistic particle on a background, use its 
stress energy tensor as a source, and solve the equations for its gravitational field. Because of 
Lorentz boosts, the field becomes completely confined to the space transverse to the direction 
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Figure 1: 2D conical bulk geometry of a thin brane vacuum. 



of motion, and hence the field equations, however formidable they may be, will linearize. A 
straightforward technique to solve them is to use the cutting and pasting technique of Dray 
and 't Hooft [65, 66], which works in our case since the sub-critical vacuum metric Eq. (4), 
with b < 1, is locally flat. So: the shocked metric will be 

ds 6 2 = Adudv - 45(u)f(x±, y)du 2 + dx]_ + dp 2 + (1 - b) 2 p 2 d<fi 2 , (5) 

where x± coordinatize the two directions along the brane transverse to the particle that 
moves along the null geodesic u — 0, and y are the coordinates on the cone. The function 
/ is the shock wave profile that we need to solve for. Then, using Eq. (5) as an ansatz, 
we substitute it into the field equations Eq. (1), with the relativistic stress energy tensor 
t^ v = -^g4uv$(u)5 ( - 2 \x±)5 fJ ' v 5 u L , of the source with 4-momentum p^ = (p, 0,0, p) included, 
and determine the equation controlling the shock wave profile /. 

Since for relativistic particles = 0, we still have Rq = jjk5^ 2 \y) and R$ = regardless 
of the source. We substitute this in the field equations along the brane worldvolume, 

MtR,\ + M\Rf v 5^ (y) = r\ 5® (y) . ( 6 ) 

The field equations in the bulk, R 6 = j^^ 2 \y), are trivially solved by Eq. (5) and so all 
the new information is completely contained in (6). The solution of these equations will be 
consistent, since one can easily verify that V^r^ = 0. Now, the only nontrivial components 
of the Ricci tensors for the metric (5) are 

R e v u = 5(u) V 4 2 /, 

RA V u = 5{u)V 2 2 f, (7) 
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where V4 2 is the full transverse Laplacian, and V2 2 its restriction on the 2D transverse plane 
on the brane, coordinatized by x±. Substituting this and the formula for t^ u in Eq. (6) we 
finally obtain the field equation for /: 

V42/ + W£ Wf ^ = W/ 2 "^ 6(2) ^ ■ (8) 

All other equations in (6) are trivially satisfied on (5). Thus we see that the problem of de- 
termining the configuration which solves (modified) gravity equations for relativistic sources, 
as before, maps onto a much simpler problem of solving a (modified) Poisson equation for 
a static charge [65]-[66],[16]. This illustrates the power of shock therapy as a diagnostic 
method. Once we have the solution of Eq. (8), we can compare it to the Aichelburg-Sexl 
solution [57] to see if, and how, the theory mimics AD General Relativity. 

Let us now solve Eq. (8), at least formally. Fourier transforming in the brane transverse 
space 

5 2 (x ± ) = j^y 2 f d 2 ke^- , / = JL J tf% Mv) e &3 ± , (9) 

using V4 2 = Vy> 2 + V2 2 , where the operator Vj/ 2 is the Laplacian on the 2D cone, and 
replacing S^(y) = 2 -K{i-b)p ^P) using the axial symmetry of (5) and the null source t^ u , 
after a little algebra we obtain 

^- k ^ = Mrhm (2p+M!k2vt)S -f- (10) 

which is very similar to the field equation of a bulk scalar with kinetic terms induced on the 
brane [15]. Note that the bulk coupling Mg~ 4 has been replaced by the effective coupling 
[(1 — 6)Mg] _1 which is stronger on the cone because of its deficit angle, representing the 
'lightning rod' amplification effect pointed out in [50]. At any rate, we can now use the same 
tricks to solve it as [15]. Ignoring the angular variable by axial symmetry, and substituting 
Whip) — D{k, p)F{k) where D obeys the Helmholtz equation 

( V 2 k 2 )D _ 2P S( P ) 

we then determine F such that the ansatz (p k = DF correctly solves Eq. (9). This yields 
F = (l — ^k 2 D(k, 0)) , and therefore the momentum space field (p k is 

Vk = gfe£> . (12) 

The Helmholtz equation is defined on the full 2D cone, with the source residing exactly on 
its vertex, and so we can write its solution D in terms of the modified Bessel function K , 
which vanishes at infinity and hence describes properly gravity localized to the brane as 

^^-W^Ml^- <13) 
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where we have explicitly used the fact that the other modified Bessel function I is equal 
to unity on the brane at p = 0, while it diverges at infinity. We point it out since it will 
arise naturally later on, although here it plays no role yet. Then, substituting this in tp^ 
and Fourier transforming it back to the transverse configuration space along the brane, and 
integrating over the angular bulk variable in the usual way 1 we finally find the solution for 
the shock wave profile on a thin 3-brane: 



This solution is very interesting - because it is pathological. Its pathologies stem from 
the fact that K diverges in the core of the brane. Indeed, K (e) ~ — ln(e) — > oo as e — > 0. 
However, along the brane, at p — 0, the divergence precisely cancels between the numerator 
and the denominator in Eq. (14). In this case, we find 



which is exactly the 4D Aichelburg-Sexl solution - at all distances along the 3-brane! This 
result would be incredibly interesting, because it occurs despite the infinite bulk space sur- 
rounding it, that could be explored by the graviton multiplet. This looks like a very efficient 
mechanism to hide extra dimensions from the gravitational exploration. However, if we 
move off the brane by an even infinitesimal amount, to p = e ^ 0, from (14) we find that the 
divergence of the denominator forces the shock wave profile immediately to zero. Thus the 
confinement of gravity to the brane is perfect - but it happens because the Bessel function 
diverges in the core, and is clearly a UV sensitive answer, as indeed noted in [15, 29]. If we 
consider a physically more realistic brane of finite thickness, the confinement may be far less 
than perfect, leading to the reopening of extra dimensions to gravity at some finite distance 
from the source [29] . Therefore, before we draw any conclusions about brane induced gravity 
in codimension> 2, we first must regulate the short distance singularities inside the brane. 

Our shock wave example points us to a very simple way how to regulate. The discontinu- 
ity of the solution is due to the divergence in the function K (p) at p = 0. Mathematically, 
the divergence of K on the source is the familiar pathology from an ill-posed problem for 
an elliptic Helmholtz equation for a pointlike source, which does not have a regular exte- 
rior solution in 2D. A well known prescription for regulating this problem is to replace the 
pointlike source by a ring-like one. We can then split the problem into a combination of 
interior and exterior problems, which have well behaved solutions for ring sources of a finite 
size. Since our pointlike source is really a particle on a codimension-2 brane, a natural way 
to replace it by a ring-like source is to consider a loop of string moving on a cylinder. If we 
then want to derive the background cylinder from a covariant action, the easiest way to do 
it is to take a 4-brane and wrap it around a circle of finite size. Such tricks were used earlier 
in [43, 51] and we employ it next. 

1 We use / = / d 2 k(fik e lk ' x± and integrate over <f) to get / = -j^ / °° dk k ipk Jo{k\x±\), as in any 
cylindrically symmetric problem in potential theory. 




(14) 




(15) 
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3 Wrapped Branes and Truncated Cones 



3.1 Making Vacua 

We wish to replace the thin brane vacuum (4) by one where the brane has a core of finite 
thickness, but its metric remains the same as (4) outside of the brane core. In this case 
the thick brane vacuum will retain some of the attractive features of the thin 3-brane, 
namely 4D flatness and the presence of 4D General Relativity, for the sub-critical case. 
Furthermore, we will find a simple model which will enable us to see what goes on as the 
tension approaches and exceeds the critical limit. Now, the reason why the metric Eq. (4) 
solved the field equations (1) was that the stress energy tensor of a tensional 3-brane has the 
form T^ v = —XS^(y) along the brane and vanished away from it. In component notation, 
this looks like T A B = -\6®(y) diag(l, 1, 1, 1, 0, 0). A very simple model of a 3-brane with a 
finite core is a cylinder: a 4-brane wrapped on a circle [43, 51]. However, when the 4-brane 
has nonzero tension, with the vacuum action S va , cunm — — J dPxy/gHXs, its stress energy tensor 
will be T A b = —X^5{p — ro) diag(l, 1,1,1, 1, 0), where we wound it around the circle of radius 
r , with p along the normal to the brane. Thus, to wrap the 4-brane into a cylinder, we 
must cancel the pressure oc A 5 in the compact direction. 

A simple way to do it is to put an axion-like field E on the 4-brane, with a vacuum action 

S_ = - J d 5 x^g- 5 (X 5 + ^djld^) , (16) 

and look for solutions where S breaks the translational invariance in the compact direction, 
but so that its stress energy tensor remains translationally symmetric. This is precisely the 
method of Scherk and Schwarz [53] for dimensionally reducing supergravities to get massive 
lower dimensional theories. The stress energy tensor obtained by varying Eq. (16) along the 
4-brane is 

T\ = -X 5 5\ + d a mz - ^5 a b g cd d c Zd d Z , (17) 

where lower case latin indices a,b, . . . run over 4-brane worldvolume values 0, 1, . . . , 4. Taking 
the fourth coordinate to be the angle on the compact circle and substituting the Scherk- 
Schwarz ansatz £ = qcp we choose q to precisely cancel T 9 ^. This requires 

A 5 = ^V , (18) 

where = \ is the inverse radius squared of the compact dimensions. Then the remaining 
components of the 4-brane stress energy become precisely T^ v = —2X 5 5^ U . Thus the brane 
source now reads T A B = —2X 5 5(p — r ) diag(l, 1, 1, 1, 0, 0). The tensor structure is precisely 
the same as in the stress energy tensor of a thin 3-brane. Now, since the thin brane is axially 
symmetric from the bulk viewpoint, we can rewrite the scalar 5-function in polar coordinates 
as S^(y) = 2~<Kp)- By inspection, we see that if we shift the argument of the remaining 

radial 5-function to p — r , it becomes S^ ick (y) = -^-^{p — r o), and so to model the thick 

f<y\ 

brane by a wrapped 4-brane, we need to identify A5 th - ck (yJ with 2X^5{p — r ). This yields 

A = 47rr A 5 , (19) 
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relating the 5D and the effective AD vacuum energies. The effective AD vacuum energy A 
contains the contributions from the classical axion field E, doubling up its value, because 
of the cancellation condition of Eq. (18). This is equivalent to the effective vacuum energy 
contributions in the string landscape [3], which is also comprised of the contributions from 
quantum fields (here, A 5 ) and classical fluxes (here oc q 2 ). In the landscape, these contri- 
butions need to be finely adjusted to yield a AD Minkowski vacuum. Here, the situation 
is similar: because this stress energy is completely localized on a codimension-1 cylindrical 
surface of finite radius, one expects that the solution representing this configuration should 
be obtained from patching locally flat 6D Minkowski spaces inside and outside of the cylin- 
der. The exterior should have a deficit angle by the construction of the thick brane stress 
energy, which shows that our 4-brane source looks precisely like the thick 3-brane. 

Indeed, the field equations for a 4-brane in a 6D bulk which include brane localized gravity 
terms are [15, 29], in Gaussian-normal gauge, and with the 4-brane residing at p = r , 

M*G 6 A B + M 3 5 G 5 a b 5 A a 5 h B 5(p - r ) = T\5 A a 5 h B 5(p - r ) , (20) 

which we apply to our cylindrical brane vacuum, imposing that the stress energy tensor is 
covariantly conserved, that requires d a d a T, = 0. It is easy to check that by axial symmetry the 
ansatz E = q(f> trivially solves the latter equation. Then, the tensor structure of the source, 
T A B = — 2\ 5 5(p — r ) diag(l, 1,1,1, 0, 0), guarantees that the flat AD metric is a solution, if 
the tension is again off-loaded into the bulk, like in the thin 3-brane case. Tracing the field 
equations (20), and using (18) we find that the condition for this is that the metric in the 
remaining two dimensions, coordinatized by the angular direction along the brane and the 
radial distance away from it has curvature 

R 2 = ^5(p - r ) . (21) 

Thus, by the arguments identical to those we invoked in the construction of the thin 3- 
brane solutions, we can choose the coordinates where the 2D metric is ds 2 = e~®dz 2 , where 
however the brane now resides at ro = \zq\ around the origin. Substituting in (21), we find 
R 2 = e~®Vy). Comparing to (21) and using the fact that V|[0(r— ro) ln(r/r )] = 5(r—r )/r, 
we infer that 

■d = 269(r - r ) ln(— ) (22) 
ro 

where 6(x) is the Heaviside step function, Q(x) = for x < and Q(x) = 1 for x > 0, and 
b is given precisely by combining Eqs. (2) and (19), 



2A 5 r 
M« 4 



(23) 



The solution is given by ds§ 2 = r]^ v dx^dx v + e 260(r r o) ln<r / r o)(^r 2 + r 2 d(p 2 ), in polar coordi- 
nates. After defining the new radial coordinate p by dp/dr = (r /r) be<r ~ r °), such that 

r fr_o\ b@ ( r ~ r °) _ &r e(r-r ) 

P l-69(r-r ) \r) 1 - bQ(r - r ) ' 1 ' 
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and 0(r — tq) = Q(p — Tq), we can rewrite the metric as 



dse 2 = ri liv dx il dx v + dp 2 + 



[1 - bQ{p - r ))p + br Q(p - r ) 



(25) 



For p < To the 2D part of the metric is ds^ 



dp 2 + p 2 c 



i.e. a flat disk, while for p > tq 



it becomes ds 2 2 = dp 2 + ((1 — b)p + br ) 2 d(f) 2 = dp 2 + (1 — b) 2 (p + ^ r ) 2 d(j) 2 , precisely the 
metric on a cone, which we see by comparing it to the thin 3-brane metric from the previous 
section. Thus the combined solution represents a truncated cone, with a flat mesa of radius 
r on top, as depicted in Fig. (2). The geometry of the brane is 4D Minkowski x circle. 




brane 



Figure 2: 2D bulk geometry of the resolved brane. 



Looking at Eq. (25) for p > r we see that the critical string tension, where the deficit 
angle becomes 2tt, is again given by b = 1, which now corresponds to the tension 

^5cr = 77-^" • (26) 
2r 

Comparing to the relation between the 5D tension A5 and the effective 4D tension A, given 
by Eq. (19) we find that the critical value of the effective 4D tension is, not surprisingly, 
X cr = 2-kMq. We can now see what happens as A — > X cr and beyond. Since we are holding 
the brane core fixed, using its own internal dynamics, all we need to do is analyze the region 
p > r of Eq. (25). When 6=1, the metric outside of the brane changes to 

ds 6 2 = rj^dx^dx" + dp 2 + r^dc/P , (27) 

which describes a cylinder M 4 x R x S 1 , where the radius of 5* 1 is the same as the brane 
radius. Now, when b = 1 — e, it's easy to see that the bulk cone looks like a sliver. In this 
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limit, the exterior bulk metric is approximately ds§ 2 = •q^dx^dx 1 ' + dp 2 + (e(p—r ) + r ) dcf 2 . 
This looks like the metric on a cylinder for r < p r /e, because the variation of the radius 
of the sliver, due to the change of radial distance p, is very small compared to the radius of 
the brane r . Near-critical branes reside inside deep bulk throats that only asymptotically 
open into conical geometries [59], see Fig (3). 




Figure 3: 2D bulk geometry of the near-critical resolved 3-brane. 



For supercritical branes, b > 1, the exterior metric is just 

ds 6 2 = r ]llu dx^dx u + dp 2 + (b- l) 2 (p - r ) 2 # 2 , (28) 

which describes a bulk geometry that near p = tq looks like a flat disk. However now the 
metric has a singularity at p = j^r > r , where the geometry tapers off to a conical 
spike. Thus when b > 1 the bulk spontaneously compactifies on a 2D teardrop, precisely 
of the shape discussed by [55]. Hence static flat branes in an infinite bulk will only occur 
for sub-critical values of the tension. In light of the discussions of [54], one expects that 
the supercritical brane solutions without bulk singularities exist, but that they describe 
nonstationary backgrounds. Examples describing topologically inflating domain walls on 
codimension-1 brane induced gravity models have been discussed recently in [56]. 

Note, that our construction is completely covariant from the bulk point of view, since 
we can define the regulator by a variational principle, from a well defined 4-brane action. 
Therefore we can explicitly check the properties of brane localized gravity, and calculate the 
crossover scale inferred in [29] exactly. 



3.2 Stress Energy Sources on Thick Branes 

Having constructed the thick brane vacua, we can turn to the description of matter sources 
that reside on branes. For simplicity we take the matter to be localized on the 4-brane, and 
described by a matter sector action 

^matter = — J d^X^fg^ £ m atter > (29) 
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where the matter couples to the induced metric on the 4-brane g^ a b- Since one of the 
dimensions along the brane is a circle of radius r , and we are interested in the behavior of 
the theory at length scales r ^> r , we dimensionally reduce the brane-localized sector on 
this direction, representing all brane fields by the standard Fourier expansion 

oo 

<M^Vo0) = E $A/>(^)e^, (30) 

n=— oo 

where M represents the quantum numbers of any particular representation on the brane. The 
Fourier coefficients then give rise to the 4D KK towers of fields, with masses M 2 = m 2 +n 2 /r 2 , 
where m are the explicit mass terms from (29), and the mass gap is 1/ro- This expansion 
also applies to the brane induced metric and its curvature. Thus, at distances r > r we can 
always model the angular dependence of any configuration as a small, Yukawa-suppressed 
correction from a KK momentum mode stretching along the compact circle. 

When we consider the classical limit of the theory, to define the gravitational sources 
that describe lumps of matter living on the brane, we can model them as loops of string 
along the cylinder. At distances r ^> r these will look just like point particles on a thin 
3-brane. To determine the stress energy tensors of such objects, that source their gravity, 
we replace (29) by a Nambu-Goto action for a string, 

Smattcr = ~H j &<Ty/l , (31) 

where \i is the mass per unit length of the string, a a are its two worldsheet coordinates, which 
we will gauge fix to be t and r o 0, and 7 Q/3 is the induced worldsheet metric. We then vary this 
action with respect to the metric of the target space on which the string moves, in this case 
the metric of the background thick brane. We then use the standard definition of the stress 
energy tensor r ab , given by the functional derivative equation 5S maUeT — \ J d 5 x ^/g5Sg a bT ab , 
to read off the stress energy tensor of the loop. Following the standard tricks for pointlike 
sources, we first rewrite (31) as 

S m = -fiJ d 2 ad 5 x v ^5 (5 \x M - x M (a)) , (32) 

where the 5-function puts the string on its shell (i.e. enforces its equation of motion). Using 

7a/3 = 9ab§£§ji and 5^f = \^lT p 5^ al3l we find 

^ = -^^^ e -^"^- (33) 

Similar variation with respect to the target space coordinates of the string leads to the string 
equation of motion, which is V a (<yuV a x b ) = 0, and is in fact completely equivalent to the 
string stress energy conservation, V a r ab = 0. 

With the gauge choice a = x° and a 1 = r o 0, and using the vacuum solution (25) as the 
target space, we get 7 a/ g = i] a p, and hence the loop of string stress energy tensor becomes 

r- = -„/ .ftrfOV - X ^)) V ^^. (34) 
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While this expression is only valid on the locally flat target (25), as we will only be interested 
in the linearized gravitational fields (which however are exact for relativistic particles!) this is 
all we shall need. If we consider static sources, the string representing only the lightest modes 
must be translationally invariant in the compact direction, because by the discussion above 
any inhomogeneities along the string can be viewed as heavy KK states at distances r » r . 
For such axially symmetric sources, only two components of (34) survive: r 00 = p5^(x), 
and r 9 ^ = — f^5^(x), which we can write more compactly as 



T a b = -n ^(3) (£) diag(l, 0,0,0,1). (35) 

The parameter p is the total rest mass per unit length of the configuration. To get the stress 
energy tensor of a relativistic string, that moves along the cylindrical brane at the speed of 
light, and has vanishing rest mass, we can simply boost (35) along the direction of motion 
x\\, go to the light cone coordinates x\\ = v + u,t = v — u and take the limit of infinite boost 
parameter while sending p — > [65]. Since the momentum of the string p is the length of 
the string 2nr multiplied by the finite limit of p cosh (3 as f3 diverges, for sources composed 
of the lightest modes as in (35) this yields, thanks to the properties of the vacuum (25), 



- 2P ^ 9au V 5(u)5^(x ± )5^J\ , t% = 0, (36) 



which is identical in form to the relativistic stress energy on thin branes, used in the previous 
section. 

To find the gravitational fields of these sources, we need to solve the brane induced 
gravity field equations (20), in Gaussian- normal gauge with the brane at y = p — r = 
and with T a b, which is given by the sum of the vacuum contribution in Eq. (17) and r a 5 
as given in (35) or (36), on the right-hand side. It will turn out to be convenient to rewrite 
these equations as the conditions on the Ricci tensor, which after tracing over the indices 
and rearranging the scalar contributions, become 



M*R 6 A B 



(t\ - Ml(R 5 a b - ^5 a b R 5 )y A a 5 b B - l -5 A B (T\ + ^Mfi? 5 



S( P -r ). (37) 



We first consider the relativistic sources obeying (36) and find the generalization of the shock 
wave solution (5), (14) on thick branes. 



4 Shocking Crossover Physics 
4.1 Setting Up the Shock Wave 

As before, we determine the shock wave solution by introducing a discontinuity in the metric 
of Eq. (25) according to the prescription of [65, 16], and demanding that the wave profile 
solves Eq. (37) with (17) and (36) as sources. We are working with sub-critical branes, 
b < 1, which admit static flat vacua. We can still, however, consider the limit b — > 1. Using 
the null coordinates u,v we require that the source moves along u = 0, and that the wave 
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profile explores all the transverse dimensions available to gravity. By axial symmetry of 
the source, however, we only need to look at the profiles f(x±, p)8(u) that depend on the 
transverse coordinates along the brane x± and the transverse distance p orthogonal to the 
brane. Thus we seek the shock wave metric in the form 



ds 6 2 = Adudv - 45(u)fdu 2 + dx{ + dp 2 + 



(l_ & 0( p _ ro )) p + &ro 0( p _ ro ) 



2 

d(t> 2 . (38) 



Now, we substitute the sum of Eqs. (17) and (36) along with the condition (18) in Eq. (37) 
and trace over the indices. We note that the stress energy of a relativistic particle (36) is 
traceless, reflecting that the source is scale-free. We can then easily check that by the form 
of the shock wave metric (38) the trace of the 6D Ricci tensor picks up contributions only 
from the p-(f> part of the metric, R e = R 2 , leading to an equation identical to Eq. (21). Thus 
as before this equation is trivially satisfied by Eq. (38), with the deficit angle parameter 
defined in (23), and so are all the tensor equations on the truncated cone p-0. 

Among the remaining tensor equations of (37), the only one which contains nontrivial 
information, as before, is oc R$" u . Indeed, Eqns. (37) on the background (38) reduce to 

MtR % \ + M\R^ V 5(p - r„) = r"„ 5(p - r ) , (39) 

and the only components of the Ricci tensors which depends on the shock wave profile are 

R 6 \ = 5(u)V 4 2 f, 

R 5 V u = S(u)V 2 2 f, (40) 

where we are explicitly using axial symmetry in the latter term when we write the Laplacian 
in only the two transverse spatial coordinates x±. Substituting this and (36) in Eq. (39), 
after a straightforward algebra we finally get 

Wf + j%WfS(p - n.) = j^^J Hp - oO ■ («) 

This equation is almost identical to Eq. (8) governing the shock wave profile on a thin 
3-brane. The main difference is that the 5-function is shifted from the origin and smeared 
over a circle p = ro, and that the Laplacian V4 2 is defined on the truncated cone (25), such 
that for axially symmetric shocks, V 4 2 / = V^ 2 / + /" + {l _ m{ pI^ P p ~l^ e(p _ n)) f\ where the 
prime denotes the p-derivative. The latter term contains a jump across the 4-brane. 

4.2 Solving for the Shock Wave Profile 

To solve (41), we again Fourier transform in the brane transverse space, using Eq. (9). This 
yields 

( V »-*>' = (2^W + I^)*- r °)- < 42 > 

In contrast to Eq. (8) here we have a single 5-function. Of course, this is because this time 
the brane is codimension-1. Hence we can use the standard methods for finding bulk field 
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wave functions in the codimension-1 brane setups: we solve Eq. (42) inside (— ) and outside 
(+) the brane, and then use the boundary conditions, which are the continuity of the 
jump of tp' k as prescribed by Gaussian pillbox integration of (42) around the brane, and the 
regularity of the solution in the center of (25) and at infinity: 

Vfc(O) = o, 
y^(oo) = 0, 



tif'(ro) - ^'(ro) = U-^ + ^V,(ro) ■ (43) 



,27rr M 6 4 M 6 4 
Away from the brane the differential equation (42) gives 

W + br ]e {p -r Q ) ^' - fcVfc = , (44) 

P + l-be(p-ro) 

with Bessel functions I (k(p+ rzgg^ )), K (k{p + ^fc"^ )) as solutions. The first two 
of the boundary conditions in Eq. (43) then pick ip k ~ I , ip k ~ K , and the third condition 
sets ip k = AI (kp)K (k{^) and = AI (kr )K (k(p + f^)). Finally the last of (43) 
fixes the coefficient A, yielding 

-4 = , , * , : r. (45) 



2TTrr 



Mlk^hikr^K^k^) + h(kr )K (k^- b )) + MlkH (kr )K (k^- b ) 



where we have used the recursion formulas for modified Bessel functions 8 z Iq(z) = I\{z) and 
d z K (z) = —Ki(z). Substituting this and the form of the solutions <pk back in the Fourier 
transform equations (9), and again using the axial symmetry in the x± plane to integrate 
over the transverse spatial angle, yielding / = ^ dk k(p k J (k\x±\), we finally obtain the 
expression for the shock wave profile on the truncated cone: 

■ 27r 2 r io Mi^k^K^k^ + hik^Koikf^)) + MikI (kr )K (kf^y 

(46) 

where 

{ Io(kp)K (krf\) , p<r ; 

I (k P< )K (k P> ) = { } \ \ 1 ~ bJ \ \ (47) 
\l (kr )K (k(p + $.)), p>r . 

This solution is clearly not divergent away from the brane, unlike (8), as long as tq > 0. 
Indeed, along the brane, p = r , it reduces to 

P r,, I (kr )K (k^- b )J (k\x ± \) 

ClrC ■ 



2vrV io Ml^kr^K^kf^+hikr^K^kf^)) + MikI (kr )K (kf^) ' 

(48) 

where all the individual contributions remain finite. 
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4.3 Crossover Physics 

Now we turn to the physics of the solution (46)- (48). It is subtle, controlled by the deficit 
angle of the background (25). Let us look at the field along the brane, (48). The shock 
wave profile in (48) is dominated by the modes with the momenta k ~ due to the 

oscillatory nature of Jo- The contributions of the modes with momenta k far from 
will interfere destructively. At transverse distances much larger than the size of the compact 
dimension, ^> r , we need to focus on the momenta for which kr <C 1. Thus we can 
always replace the terms oc I n (kr ) by their small argument expansion. For sub-critical 
branes, 1 — b ~ 0(1), and we can likewise replace any K n (kj\) by their small argument 
expansion too. However, in the near-critical limit the deficit angle approaches 2tt and so 
|1 — 6| C 1. Hence the argument of K (kj^) may be very large even when kr <C 1. Thus 
we must consider the near-critical branes very carefully. 

We start with sub-critical branes where kr <C 1 — b, and we can approximate all Bessel 
functions by their small argument form. Using I v {z) — > r ^ +1 ^ (§)", K v (z) — > ^^(f) I/ for 
v > and K (z) — > \n(2/z), the shock along the brane then becomes 

r> r°° In ^r— — Jo(k\x±\) 

If k ~ 1 / 1 | is large enough for the second term in the denominator to dominate, which - 
in this approximation - happens when 

(!- & ) M 6 4 



k 2 > K(k) = „ , , (50) 



Mfr In 



2(1-6) 



l 5' 

the logs in the integrand cancel, and the shock wave profile is 

/( ^' ro) "-2^Mf^A c dk ^— = ^Mi ln { kcM ) + ---- (51) 

where we have used Gauss law to normalize the lightest KK modes in the expansion of 
the brane induced curvature on the circle of radius tq to Mj = 2nMj?ro. This is - to the 
leading order - precisely the 4D Aichelburg-Sexl solution, with the log profile normalized to 
the inverse of the critical momentum k c given in (50), as opposed to some arbitrary short 
distance cutoff I. 

For sub-critical branes when 1 — b ~ 0(1), this is the end of the 4D story: all distances 
longer than r fall in the kr 1 regime, and the gravitational field looks 4D, like (51), only 
when r < r c = l/k c . In this, generic, case, the crossover scale is qualitatively very different 

from the see-saw scale of [29]. Indeed, using Gauss law Mf = we can rewrite l/k% from 
(50) as 

r 2 (k) M l iJ 2 ^-^ (52 ) 

or, ignoring the log, r c ~ ttt 3 ^, where Mg eff = (1 — &)Mg. Thus for sub-critical branes 

6 cff 

the crossover scale is set by the 'naive' ratio of brane and effective bulk Planck scales, and 
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depends on the UV cutoff r of the brane core only through the logarithmic correction 
reminiscent of 'running' in real space, which was discussed for locally localized gravity in 
[63]. This shows that the crossover weakly depends on the momentum, and is very different 
from the see-saw mechanism of [29]. We also see that the bulk gravitational coupling 1/Mg 
is amplified by the 'lightning rod' factor (1 — b)^ 1 of [50]. This 'lightning rod' amplification 
is intuitively clear, because gravitational field lines spread more slowly on a cone than on 
a plane, at a fixed distance from a source, and hence gravity must look stronger on a cone 
[50]. It plays an even more important role in the near-critical limit, as we will see below. 

At distances larger than r c , gravity eventually changes to 6D. For sub-critical branes, 
this transition actually looks like a rapid 'running' of the coupling of individual graviton 
modes in the resonance that impersonates General Relativity. We can glimpse this, for 
example, from Eq. (49), which shows that the effective coupling controlling the transfer of 
the momentum k is 

■2(1-6)- 



MU (l-b)M£ ln t kr J" (53) 

In contrast to the 5D case [16], the running here is logarithmic, because there are two extra 
dimensions as probed by the localized masses on the wrapped brane. To recover the full 6D 
form of the wave, we can take (47) or (48), and taking the limit of very large distance, ignore 
the second term in the denominators of these integrals. In this case, 

f-- 2 n*(l-b)Mi f dkkI o( kr o)K (k(p+- [ ^))J (k\x ± \) + ..., (54) 

where the integral can be done in closed form, leading to [50] 

f ~ P \ |_ (KK) 

J 2n 2 {l-b)M$\x ± \ 2 + p 2 1 J 

This is the 6D gravitational shock wave, constructed in [58], with the 'lightning rod' term 
accounting for the brane tension as found in [50]. 

Let us now consider the near-critical limit, b — > 1, which requires particular care. In this 
limit, Eq. (52) suggests that the crossover scale diverges, indicating that gravity remains AD 
out to extremely large distances, regardless of the brane thickness. However our discussion 
of the near-critical limit in Sec. (3), in particular interpretation of Eq. (27), has already 
taught us otherwise. As we have seen there, in the near-critical limit the bulk around the 
brane compactifies to a long cylinder, which opens up into a cone only very far from the 
brane. Bulk gravity in the throat is in the 5D regime, which separates the 4D and 6D ones. 
Thus the crossover scale beyond which gravity is not 4D is really the demarkation between 
the 4D regime, and 5D gravity which lives inside the throat. Looking at the shock wave 
solution (48), we see that while the small argument expansion for I n (z) always suffices for 
large distances from the source, in the near-critical limit we must treat the functions K n {z) 
carefully for the momenta in the regime 1 — b <C kr <C I. Indeed, this corresponds to length 
scales £ < ^ which are shorter than the length of the conical throat. So setting I (kr ) — > 1 
and ii(Axo) — > in (48), we notice that where the terms oc Mf dominate, gravity will remain 
4D because the factors of K (kj^) exactly cancel in the integrand. This happens for the 
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M 4 Ki(k r " ) 

momenta k which satisfy k > -rA ,, % b . . In the near-critical limit this can occur when 

the argument of Bessel functions is much larger than unity, so that the ratio Ki/K will 
be of the order of unity! Indeed, as we decrease kr , going out to larger distances, K 1 will 
remain close to K Q : since K\(z) = —d z K (z) and now, since b — > 1, we have large values of 
the argument for which K (z) — > y^e -2 , so that fTi — > fTo- The transition can occur only 
once, because K n are monotonic, and for our approximation to be valid it must happen when 
Y-b — ^(1) " otherwise the situation will be essentially the same as in a generic sub-critical 

M 4 

case. In that case the crossover scale saturates to k c = jffi, or therefore 



Ml 



(56) 



Clearly, this is valid as long as r c < -f^. The scale r c is exactly the see-saw scale of [29]: 
again using Gauss's law for the AD Planck mass, Mf = 27rMfr , we find 

r - — (57) 

This demystifies the appearance of the see-saw effect, showing that at least with our regu- 
larization, it does happen naturally, but only for near-critical branes whose deficit angle is 
very close to 2n, so that they still are flat and static. 

As an explicit check that inside the conical throat gravity looks 5D, we can take the limit 
of our solution for a fixed distance and formally take the limit b — > 1 while holding tq 
fixed. Then if we take the solution in Eq. (48), in the limit b — > 1 we can replace K" n 's by 
their large argument expansion for any finite momentum. After cancelling the like terms, 
and taking the thin brane limit, we find, using footnote (2), 

fr x P [°° Jo(fc|£i|) _ 9 f°° _fk ef^ 

J{x±,P) 27r2 J Q aK M ± rQ + M 3 rok ^y_ oo ( 27r )2 27cM£r k + 2nMlr k 2 ' 1 ' 

This is precisely the shock wave profile on a flat thin 3-brane in 5D, on the normal branch, 
with brane and bulk Planck masses M| e // = 27rM|r and M^ e ^ = 27rM|r , respectively. It 
can be obtained as a limit of the normal branch solutions of [16] on a tensionless brane. 

So to sum up, we find that for distances r < r < r c the shock wave reduces to the 
Aichelburg-Sexl solution [57], and becomes its higher-dimensional generalization when r > r c . 
The crossover scale for sub-critical and near-critical branes, beyond which gravity is not 4D 
anymore, is given by 



M 



2 _ I 2n(l-b)Ml 



In 



2(l-fc 



1 U 1 1 

4w 2 Mgri ' — > 1 



6<1 ; 

(59) 



J 6 ' 

only in the latter case agreeing with the see-saw scale discussed in [29]. The origin of the 
see-saw mechanism is now clear: in the near-critical limit, the brane resides in a very deep 
throat generated by its own tension. Thus its localized 4D gravity first changes to 5D gravity 
residing inside the throat, and only later to full 6D gravity [59], at distances p 3> r /(l — b) 
(see discussion of Eq. (27)). This shows how the see-saw mechanism emerges. 
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5 Linearized Gravity 



5.1 Perturbations and Their Gauge Symmetries 

We can solve exactly the field equations (1) for localized matter sources only in the relativistic 
limit. For nonrelativistic sources, we can try to solve them perturbatively, assuming that 
the solutions converge to the background vacuum far from the sources. So starting with a 
background ds^ = qab dx A dx B metric, we define the metric perturbations Kab as is usual, 
by the substitution qab — > Qab + ^ab- In addition, we must also perturb the brane location 
r o —> r o + £ and any nontrivial field configurations on the brane, in our case the axion E — > 
E + a. Before substituting these expansions in the field equations (1), we must account for 
the gauge symmetries of the perturbations, to identify physically meaningful variables which 
cannot be undone with diffeomorphisms. So suppose that we transform x A — > x A = x A + x A - 
A straightforward calculation then shows that the perturbations transform according to 

h A B ->h~AB = hAB-^AXB-^BXA, £ -> £ = £ + X"\p=ro , O -> O = a - X A \ P =r V^S , 

(60) 

where the Va and raising and lowering of indices are defined with respect to the background 
metric qab- On the static thick brane background (25), we can immediately compute the 
transformations (60), splitting them up with respect to different representations of the 4D 
Lorentz symmetry. They are 





hpu dpXu 




h<f,<p 


= h^ — 2d<px<p ' 


-«y, 




= hp^ - dpX<t> - 




h^p 


= hpp dpX P ~~ 


Xfi i 


hpcj, 


= h p4> - d^x 9 ~ 


X4>' + X<t> 


hpp 


= hpp 1x P i 






= e + x p (r ), 




a 


= ^-gxVo), 





where 

a = g H = ((1 - 60(p - r ))p + 6r 6(p - r )) 2 , (62) 

is the area of the disk at fixed p in Eq. (25), and x p ( r o) = X p \p=r - We need to carefully 
raise and lower the angular index 0, since although the brane is intrinsically flat, the full 
geometry isn't. Note, that while a is continuous, ^= = 2(1 — bQ(p — ro)), so that on the 
brane, the discontinuity, defined for any function as the jump of its value across the brane, 
Af = f\ + - f\_ = f(r + ) - f(r -), is = -4^, where we have used Eq. (23). 

Using these transformations we can always gauge fix the metric perturbed by axially 
symmetric sources to brane-fixed Gaussian-normal gauge, in the linear order. The key is to 
be careful about what axial symmetry means: the vector field must be a Killing vector 
of the metric, unperturbed as well as perturbed. Hence any tensor, including the metric, 
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must be Lie-derived by to zero. Now, using (61), we can readily pick x p an d X<t> to set 
h pp and h pi p to zero, respectively, treating the corresponding transformation rules in (61) as 
differential equations defining x p an d X</>- We can always choose the boundary condition 
X p ( r o) — 0, opting to keep the perturbation of the brane location £, or brane bending, to 
remain unchanged, for the moment. Next, we can pick a Xy. — X/u(p, 0) t° se t ^Vp an d 
to zero too. This may seem to induce angular dependence of by the first of Eqs. (61). 
However, it's easy to check that in the coordinates where h pp = h^ = = hp P = the 
Killing conditions enforcing axial symmetry imply that d$hp V = d^h^ = 0. Thus the metric 
perturbations Jiab are ^-independent. In fact, this proof extends beyond linear perturbation 
theory, as the reader can check by consulting [64]. Notice that we have now completely 
used up x*( r o) - we can't change a without also changing h p ^ as is clear from (61). This 
shows that the lightest in the tower of the KK vectors h pa eats the brane-localized axion a 
to become massive. Alternatively, when we pick the gauge h pa = 0, as we did here, we are 
separating out the scalar a as the Stiickelberg field of h pa . 

Dropping the overbars, note that in any Gaussian-normal gauge h pp = h p $ = = 
hpp = 0, we still have the residual gauge freedom generated by Xn-, which are independent of 
p and 0, and by X p ( r o)- These transformations do not change our gauge. They do change 
hp U , h^, £. Transformations generated by x p are the standard 4D diffeomorphisms, and we 
will find their uses later on. Transformations generated by x p ( r o) allow us to remove the 
brane bending, nailing the brane straight to its background value, p = r . In this gauge 
we can quickly find the field equations for the perturbations, by directly substituting the 
gauge-fixed perturbed metric into the field equations (1) or (37). 

Going to this gauge requires a little care, however. As we have noted above, — is 
discontinuous on the brane in the background (25). Hence when x p { r o) 7^ 0, the gauge 
transformation (61) changes the metric perturbation h </,</, dis continuously. On the other 
hand, in brane-fixed Gaussian-normal gauge the full field configuration has discontinuities 
on the straight brane, proportional to its the stress-energy. These must remain encoded in the 
5-function sources in (37). Hence the metric perturbations in brane-fixed Gaussian-normal 
gauge must be continuous across the brane, and their derivatives may have at most a finite 
jump there. In any other gauge found by using a transformation (61) with x p ( r o) 0> the 
metric perturbation h </,</, will appear discontinuous, so that the brane matching conditions on 
the perturbations can only be found by carefully implementing the Israel junction conditions 
(for an illustration, see e.g. [51]). So to get the field equations, we choose to work in brane- 
fixed Gaussian-normal gauge, where we get the brane matching conditions by using Gaussian 
pillbox integration around the brane, and relating the normal derivative discontinuities to 
5-sources in the field equations. This discussion applies completely generally to any model 
where a codimension-2 brane is resolved by a cylindrical codimension-1 brane. 

Nevertheless, it turns out the bulk field equations for perturbations are easier to solve 
by stepping out of the brane-fixed Gaussian-normal gauge, to a gauge where the p-<p cone 
appears conformally flat, and explicitly keeping the brane bending £. The latter is similar to 
what happens in other codimension-1 cases [67], except that here the brane bending must be 
turned on even in the 4D vacuum (25). The reason is that the 4D vacuum backgrounds (25) 
are not codimension-1 vacua: the axion flux breaks the symmetry between the longitudinal 
brane directions. Brane bending probes this 'vacuum' anisotropy. Using the conformally 
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flat p-(j) metric, on the other hand, simplifies the bulk field equations. For later use, we list 
here how to relate such gauges. In a Gaussian-normal gauge the only nonvanishing metric 
perturbation is h^, now denoted with the caret. Defining $ by = «$, the perturbed 
metric is 

ds 6 2 = (77^ + h^) dx^dx u + dp 2 + (1 + <f>) a(p) d(j) 2 , (63) 

with a defined in (62). Let the brane bending term £ in this gauge be non-zero. To go to 
a gauge where it is zero, we apply (61), choosing x p ( r o) = ~C That takes us to the gauge 
£ = 0, where the metric has the same form as (63), but with 

$ = 1. + ^, (64) 
a 

instead of l>. This shows how $ 'ate' £ in brane-fixed Gaussian-normal gauge. On the 
other hand, to relate a general Gaussian-normal gauge to a conformal gauge, we use the 
relationship 

a' ff 

$ = $ - — / dp$, (65) 
2a Jr 

which is a gauge transformation of the form (61) with x p = \ It dp&, and other gauge 
parameters set to zero. For this gauge transformation, x p ( r o) = 0) an d therefore it does not 
change the discontinuities of the perturbations in different gauges. Now, to see which gauge 
this transformation takes us from, we need to invert this integral equation, and determine 
$. That is straightforward. By definition of x p ■, we have $ = \x p \ so we can view (65) 
as a differential equation for x p , with the boundary condition x p ( r o) = 0. The equation is 
X p ' — f^X p — an d the solution satisfying the boundary condition is x p — V® fr dp-^- 
Thus 

(66) 

Then (65) and (61) imply that in the original gauge = y/a$. Furthermore, from h pp = 
and Eqs. (65) and $ = \x p ' we see that h pp = ^/a§ too. So indeed, the gauge transformation 
(61) defined by (65), (66) relates the general Gaussian-normal perturbed metric (63) to the 
metric 

ds & 2 = (rj^ + V) dafdx" + (1 + $) (dp 2 + a(p) d<p 2 ) , (67) 

If in addition we want to relate brane-fixed Gaussian-normal gauge perturbation to the 2D 
conformal gauge perturbation (67), we can combine together the two cases. Then, (64) and 
(65) yield 

a' fP a' 
$ = $-—/ dp^ + -i ) (68) 
la. Jr a 

relating the metric perturbation $ in a 2D conformal gauge with brane bending £ to the 
metric perturbation $ in brane-fixed Gaussian-normal gauge. Since $ must be continuous 
on the brane, we see from (68) that <3> is not. Instead, it jumps across the brane by A$ = 
— A(^-)£ = 4-^r£. These formulae will come in handy later on. 
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5.2 Linearized Field Equations and Helicity Decomposition 

We find the field equations for perturbations with axially symmetric sources by substituting 
the perturbed metric in the brane-fixed Gaussian-normal gauge, 

ds 6 2 = (rj^ + huJ) dx»dx v + dp 2 + (1 + $) a(p) d(f) 2 , (69) 

with C, = 0, into (37), and expanding to the linear order in the perturbation $. Since the 
field equation V a V a S = reduces to the source-free 4D Klein-Gordon equation for the axion 
perturbation a: d 4 2 a = 0, in the linear order the axion is not sourced by axially symmetric 
matter distributions. We can simply set a to zero from now on. This is self-consistent at 
distances \x\,p r . In this limit we can dimensionally reduce the brane theory to a set of 
KK towers with a mass gap oc l/r . Since a background on the brane which is not axially 
symmetric must involve higher-level KK states, to encode the angular variation around the 
brane, at low energies E <C l/r we can only excite the lightest states, and the backgrounds 
will automatically be axially symmetric. In fact, since the axion a really is the Stuckelberg 
field of the vector h pa , it will only be turned on when the vector field is turned on, which 
means when there is a non-trivial current that can source h pa . 

Evaluating the curvature tensors and the stress-energy terms in (37) we find the field 
equations for perturbations in this limit. They are the momentum constraint, coming from 
the p-p components, 

dx(h\-5\(K + & + ?U)) = 0, (70) 

the trace equation, 

K + & + ^ = ^4 - ^ 4 2 l> + 3A 5 $ + ™£-{&d h' x - d A 2 h A ))5{p - r ) , (71) 



the 'radion field' equation, coming from the 0-0 component, 



a ~. a 
— $' H . 

a 2a 

1 / „ , * 3Mf „ ~ _ ~ M| 



2M? 



(t\ - Ar% - ^d 4 2 $ + 7A 5 $ - -f{d x d a h° x - d 4 2 h 4 ))5(p - r ) , (72) 



'6 

and finally the 4D tensor equations, 



d^d x h\ + dM^ - V 6 2 /^ - d"d u h 4 - d»d„§ 



2 

Ml 
2 



1 M 1 A 

t\ - -r\5^ + -^(9^$ - -<9 4 2 <M^) + -i&^u - (73) 
4 2 4 4 



-f-frdxhK + d v d x w x - d 4 2 h\ - d»d v h 4 - -5\{d x d a h\ - d A 2 h 4 )) 



5(p-r ). 



Here h 4 = ^ is the trace of the AD perturbation h pu , the sources r a b describing the brane 
matter perturbations naturally decompose as r^ v and r^, and r a a = r M M + the full 5D 
trace of T a b . The covariant d'Alembertian V§ 2 is defined as before, with respect to the full 
background metric (25). 
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This system contains redundancies, which we need to disentangle away. First, we need to 
break down the perturbation h^ u into irreducible 4D representations of the Lorentz group, 
which is the symmetry of the background (25). Using the flat limit of the decomposition 
theorem in [251, we can write 



(74) 



0, with 5 helicities, 



are two scalars. 



where 7^ is a transverse-traceless (TT) tensor satisfying d^)^ h 

is a Lorenz-gauge vector, d^A^ 1 = 0, with 3 helicities, and \I/ and 
The total number of degrees of freedom adds up to ten, which is the number of independent 
components of a symmetric 4x4 tensor. Once we recall that these modes are all dependent 
on p and, in principle, could also depend on 0, we really end up with KK towers of degrees 
of freedom classified by their 4D masses in addition to their helicities. Of course, not all 
of these modes may end up being physical; some could be pure gauge modes, and one has 
to carefully check that using the field equations and the gauge symmetries of the problem. 
As we will only work with axially symmetric perturbations, our task will be considerably 
simpler. 

The field equations (70)- (73) for the linearized theory also decompose by helicities due to 
AD Lorentz invariance of the background. A simple way to separate out the equations is to 



define the transverse-traceless tensor projection operator K^ va p, such that K^ v al3 h, 
This operator must satisfy K^ ua p 
solution is 



la/3 Kfj,i/- 

K Viiaj3 = and K^^ afS = d^K^ vaj3 = 0. The unique 



+ 2 



3/ , \ 0a9 

d u d a . d u d a d^dp dj)p 



- Vaf>- q 2 
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<9 4 4 



(75) 



Clearly, this operator annihilates the vector and scalar terms in (74), and when we apply it 
to Eq. (73), it projects out the TT contributions. Since the vacuum (25) is 4_D-flat, stress 
energy conservation implies = 0, and so 



a/3 



d»d v 
~~d~J 



T 



(76) 



Hence our final form of the TT-tensor field equation is 



M 3 

V 6 V, + ^4 2 >fu6(p - r ) = 



Mi 



~~d~J 



6{p-r ). 



(77) 



Subtracting the TT-tensor equation from (73) leads to the system which contains only the 
vector and scalar contributions. We can extract the vector in a similar way as the TT-tensor, 
by defining the vector projection operator V a p, and acting on the momentum constraint (70). 
Since the backgrounds (25) are Lorentz-invariant, in the linear order of perturbation theory 
vector sources vanish and the vector modes decouple from the matter distribution. Hence 
we can set the vector field in (74) to zero. 
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This leaves us with scalars. In brane-fixed Gaussian-normal gauge, they are <£>, * and /i 4 , 
obeying the field equations (70)-(72) and what remains of (73) after we subtract out (77). 
To pick out the minimal set of independent equations for scalars from the system (70)- (73), 
we define X = \{h± — <9 4 2 \I/), write the scalar part of the metric perturbation as 

*W V = &>d v * + Xb\ , (78) 

and substitute it into the difference of Eqs. (73) and (77). The result is a combination of two 
independent AD tensors, diagonal oc and off-diagonal oc d^d u . Hence the coefficients of 
these tensors must separately vanish. Likewise, we could write the AD trace of this equation. 
However it is a trivial linear combination of the diagonal and off-diagonal terms, and so we 
can ignore it immediately. Then, combining these equations with (70)-(72), and substituting 
(78) for h^, we get 

d x (3X' + & + ^$) =0, (79) 



$" + + 9 4 V + 4X" = + r% + 3A 5 $ - ^U 4 2 ($ + 3X))*(p - r„) , (80) 



'6 

9 4 2$ + $" + + —(d A 2 m' + ax') 



2M f 4 1 ' " 

+ |V + S + 2* = ^(^|^- 2 (r%) - M 5 3 (l» + 2X])8{p - r ) , (82) 
X" + ^X' + «9 4 2 X = (J(r% - 3t%) + 2A 5 $ - M 5 3 d 4 2 ($ - - r ) , (83) 

where we have decomposed the full d'Alembertian as V6 2 = <9 4 2 + V| and used the formula 
for axially symmetric fields V| = j^2 + f^^- We have replaced the 5-D stress-energy trace by 
T a a = t^V + t'V In Eq. (79), we can immediately omit the AD divergence d\ because the AD 
background is flat Minkowski. In Eq. (82), <9 4 ~ 2 acting on r M M can be simply understood as a 
multiplication of the Fourier components of by — 1/A: 2 , as we will do later on. Finally, we 
note that the bulk terms in Eq. (83) depend only on X, and hence can be readily integrated. 
Thus we should keep it in our final set of independent field equations. 

We find that if we add Eqs. (80) and (81), and subtract from the result <9 4 2 x Eq. (82) 
and Ax Eq. (83), the result is a linear combination of the p-derivative of Eq. (79) and Eq. 
(83). Next, we can similarly start with Eq. (80) but now subtract Eq. (81), <9 4 2 x Eq. (82) 
and Ax Eq. (83). The result is very simple: 

<9 4 2 ($ + 3X) + — (<9 4 V + AX') = . (84) 

In this equation, the sources on the right hand-side completely cancel, since they are local- 
ized on the brane. The reason is that this is the perturbed gravitational equation in the 
radial direction, which is sourceless in brane-fixed Gaussian-normal gauge. This equation is 



' r%- 3t% + 7A 5 $ - ^9 4 2 ($ - X) )S(p - r ) , 
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particularly useful when we look for the solutions in the bulk, as we will see below. Further, 
we can also check, by using the p-derivative of (79), that the bulk part of Eq. (82) reduces 
precisely to Eq. (84). 

These observations teach us that in the bulk we can trade Eqs. (80), (81) and (82) for 
one independent combination, which, with wisdom after the fact, we choose to be precisely 
Eq. (84). This, however, is slightly tricky on the brane itself. A product of functions which 
are discontinuous on the brane appears in Eq. (84): and d 4 2 W + 4X'. H ence although this 
equation contains no terms oc S(p — r ), it still gives some information about the boundary 
behavior of the scalars, albeit weaker. It requires that A(^(a 4 2 ^'+4X')) =0. Although the 
factors are bounded, because they are discontinuous this condition - by itself - is degenerate: 
it relates the jumps of the functions to their values on different sides of the brane. So using it 
alone we can't extract unambiguously the junction conditions for c^ 2 ^' + AX'. However, the 
degeneracies in the boundary information in (84) are completely artificial. They are easily 
lifted by going back to any of the boundary conditions encoded in the (^-function terms in 
either of Eqs. (80), (81) and (82). From examining these equations it is clear that the most 
useful form of the boundary condition is obtained from Eq. (82), because the boundary 
condition for $ is tied to the one for X by Eq. (79). Thus, our independent equations in 
the bulk are (79), (83) and (84), supplemented with the boundary conditions obtained by 
Gaussian pillbox integration of (82) and (83) around the brane. Further, we can immediately 
integrate (79) once, with respect to x x , and ignore the non-normalizable homogeneous mode 
on the brane. With this, finally, our system of equations reduces to the bulk system 



a' 



3X' + $' + — $ = , (85) 
2a 



d A 2 ($ + 3X) + — (<9 4 V + AX') = , (86) 



2a 



a' 



X" + —X' + d 4 2 X = , (87) 
2a 

AX' = - J^H^ - 3t%) + 2A 5 $ - M 5 W($ - X)) , (88) 
M 3 ~ 

A*' = -^($ + 2X), (89) 

where the last two encode the boundary conditions, equating the jumps of bulk derivatives 
to the brane stress-energy. This set of equations is fairly straightforward to solve: Eq. (85) 
relates $ to X, which is determined by Eq. (87) and boundary condition (88). Then Eqs. 
(86) and (89) determine the remaining variable, 



5.3 Solutions 

We now turn to finding the solutions of the system (77), (85)-(89). The most interesting 
solutions are those which describe the long range fields of static sources on the brane. They 
are our key probe of what kind of gravity the theory reduces to at long distances. One 
could quite straightforwardly also determine the spectrum of the theory around the vacuum, 
when the stress-energy sources r a b = 0. We have done it, to find that the spectrum contains 
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a continuum of massive tensors and a scalar, fully described by one independent variable 
X. In this section we will give the solutions of linearized theory which apply to a general 
case of arbitrary time dependent solutions. However what is really important is to see how 
these modes couple to matter in order to determine just how dangerous, or not, the extra 
modes may be at long distances. For this reason, we will focus on the long range static 
fields, completely determined by the static rings of matter on the brane. We will also write 
down the formal solution for the full metric perturbation /i M „, which will enable us to see 
the tensor structure of full long range gravity in the linearized theory. 

Let us start with the field equation governing the tensor modes, Eq. (77). This equation 
can be solved very easily using the standard Green's function techniques. Expanding out the 
bulk d'Alembertian V6 2 , and Fourier-transforming over the 4D variables x^, while keeping 
in mind that the fields are axially symmetric, to 7^(2:, p) = J (f^r7^(&,p) exp(ik ■ x) and 

r a b(x, p) = J -f0^r a b (k, p) exp(ik ■ x), we can rewrite the TT-tensor equation as 



7V- 
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5(p-r ) . (90) 



The solution of this equation can be written in terms of the Green's function Q TT (k,p), 
defined by the equation 



Gtt" + 



^-Q JT ' - k 2 g TT - ^jxk 2 g TT S(p - r ) 
la MS 
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m 



5(p-r ). 



(91) 



Then the solution is simply 



Y„(k,p) = -2M*,p) 



k^k 



k 2 



v \ a 
I r 



(92) 



where we postpone the detailed determination of Qyy until later. For now, it is enough to see 
that it exists, after we specify appropriate boundary conditions at bulk infinity and inside 
the wrapped brane, which ensure the regularity of the solution. 

The scalar story is more intricate. The reason is that the backgrounds (25), which we take 
as 4D vacua of the theory, are not vacuum solutions from either the bulk or the wrapped 
4-brane theory points of view, because they contain a non-trivial S field contribution to 
stress-energy on the brane. As we discussed, that is necessary to wrap the brane on a circle 
and make it look like a codimension-2 object. However, because of this, the brane bending 
term will always appear in perturbation theory around the wrapped vacuum, even if the 
matter perturbations r a b vanish, in contrast to the simpler codimension-1 cases studied in 
[67]. This is interesting because it lends to a difference in the couplings of scalar metric 
perturbations, and the helicity-0 modes from the massive graviton multiplets at the same 
mass levels. 

So let us construct the solution of (85)- (89). We first Fourier-transform all the fields over 
the 4D x^ space, as in the case of TT-tensors above. Then we take a slightly roundabout 
approach, due to the presence of the background flux of E, reflected by the tension-dependent 
terms in (85)-(89) (explicit, as well as implicit in ^). We need to carefully identify how they 
mix different scalar fields. To this end it is useful to recall the gauge transformations (68), 
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and explicitly restore the brane bending term £ which keeps the brane- localized terms in 
check. Using (68) and the background equation (sl)' = — — Aj^5(p — r ), we can check 

that $' + ^<l = Substituting this in Eq. (85), we find 3X' + & = 4^£<5(p-r ), which is 

6 

very easy to integrate over p: we get 3X + $ = 4-^£©(y) +F, where F is a Fourier-transform 

6 

of some function F(x). However, we recall that this equation was already integrated once 
over x^, which implies that d^F = for consistency, or F = const.. However, any such 
integration constant can be readily removed by a residual AD diffeomorphism of the form 
X M = ^x^ as one can immediately check. Thus we can set F = 0, to find 



2 

4A 
~M, 



$ = -3X + T7 §£©(p-r ). (93) 



Again using (68) and noting the integral identity dpQ(p — r ) = (p — r o )0(p — r ), we 
can write the solution for $: 

3a' fp . „ a' 4 A, 



2a Jr a M| 



(94) 



In particular, on the brane it is given by $(r ) = — 3X(r ) + 



Next we solve for by substituting Eq. (94) into Eq. (86), and defining 

dpX +2£-4 

With this, Eq. (86) becomes 



T' = *' + 3 f P dpX + 2^-A^^e(p-r ). (95) 

Jrn Ma 



l^'-f)-^--)- < 96 » 



These two equations are straightforward to integrate. After a simple algebra we find 

* = ¥ x ~ 3 />/> A '<"' " 2 "« + W Cwj e(p - ro) ( 2( " " " 4 ) ' (97) 



where we have subtracted, using another residual 4D diffeomorphism, an arbitrary 4D wave 
which arises as the Fourier transform of the constants of integrations that yield (93). On 
the brane, this yields \l/ = — 2r £. Before proceeding, note that both solutions for $ 
and *f>, Eqs. (94) and (97), remain indifferent to the presence of brane matter sources, which 
are completely encoded in the bulk field X and the brane bending £. Once we find solutions 
for X and £ in the presence of sources, we can just substitute them back into (94) and (97) 
and have $ and ^ as well. This greatly facilitates integration, and is the reason behind our 
choice of Eqs. (85)- (89) for the description of the scalar sector. 

So let us determine X and £. First, we solve for £ using the boundary conditions (88) 
and (89) and solutions (94) and (97). Using the expression for $ on the brane given above 
we can eliminate $ from the boundary conditions. Then we use the solution (89) to relate 
A*' to AX', which gives 

A^' = ^-AX' Tph^Z- (98) 
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Finally, using (88) to evaluate AX', we compare (89) and (98), and compute £ in terms of 
the value of X on the brane. A simple albeit tedious algebra yields 



,2 | bM* 
" ~T Mfr { 

2bM£ 1 - b + /c 2 r 2 " ' 6M 6 4 1 - 6 + £; 2 r 2 ' 



3(1 -6)Mfr 2 k +M^ t (1 - 6)rg r% 



where we have used (^)|+ = ^ — ff§, obtained from (62), and resorted to Eq. (23) to 
reintroduce b = ^pr- We notice, that in the thin brane limit, r — > 0, the brane bending 

6 

term formally vanishes, if we hold all other terms in Eq. (99) fixed. Now we can finally go 
after X . Using the expressions for $ = — 3X(r ) + on the brane, and with £ given in Eq. 
(99), we can eliminate them from the boundary condition (88) in favor of X . This yields 



M 6 4 V 4M 5 3 r v 7 6M 6 4 V ° &M 6 4 1 - b + ife 2 rg 

where 



46M 6 4 1 -6 + k 2 r 2 Q V M 5 3 r , 

is the correction to the brane kinetic term of the scalar X originating from the brane bending 
contributions. The field equation in the bulk (87) is very simple. If we rewrite it with the 
boundary condition (100) restored as a 5-function term, in the Schrodinger form, we can 
compare it to the TT-tensor field equations (90). The equation is 



i ( T%+ wiziz^ * w (p _ ro) . (102) 




6M 6 4 V bM£ l-b + k 2 r 2 Q 



This equation is very similar to the TT-tensor equation (90), with the differences being in 
the brane-localized terms. We will comment on them below. Here, we note that we can now 
exploit the formal similarity with TT-tensor equation and use the Green's function methods 
again, to find the solution. This time we define the scalar Green's function Qx(k, p) by 

ft. + - fa - g(p + ™£ -«*•)) W - r„) = ±S( P - ro) . (103) 



Formally the solution is 



x(k,p) = ~g 1 (k, P ) 



T a+ bMi 1 - b + k 2 r 2 T 



(104) 



where again we defer the detailed determination of Qx until later. 

Let us now briefly outline here the procedure for extracting the long range physical fields 
from the solutions. We recall that by axial symmetry of the configurations which we are 
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exploring, at distances i > r , we can dimensionally reduce the theory on the compact 
circle. Hence the long distance dynamics will really be described by a scalar-tensor theory, 

1 /2 

where the Brans-Dicke-like scalar field, or the radion, is . In addition, in the tensor 
sector there will be the helicity-0 mode resonance as well as the helicity-2 ones. The radion 
g 1 ^ already renders the effective gravitational strength as evaluated from the brane-induced 
gravity action field-dependent. Indeed, at the linearized level about the vacuum, we see 
from Eq. (69) that since we are working in Gaussian-normal brane-fixed gauge, the scalar 
field will be Xbd — (1 + ^)ct 1//2 (p), where a 1 ^ 2 (p) is its background value, and a l ^ 2 {p) * the 
perturbation sourced by the rings of matter on the wrapped brane. Clearly, the background 
variation of Xbd retains the memory of the compactified dimension, which opens up as 
one moves away from the cylindrical brane, as discussed recently in [59]. Nevertheless, the 
description in terms of the reduced theory is still perfectly justified at the linearized level 
with axisymmetric sources. It is then useful to go to the analogue of the effective 'Einstein' 
frame for the dimensionally reduced perturbation theory, in order to separate out graviton 
and scalar forces, and compare the linearized fields to 4D General Relativity. In doing so, 
we will ignore the prefactors a 1 / 2 as they are fully absorbed in the background solution, 
and only remove the perturbations $ from the effective Planck scales, because we are most 
interested in the effective theory along the brane. Thus the radial variation of the effective 
bulk Planck scale is not important at the linearized level of perturbation theory. So we 
proceed by noting that the graviton kinetic terms in the full 6D theory come from 

r M 4 r M 3 

S 6D E> / d 6 x^^R 6 + / d b x^g~ 5 ^R 5 . (105) 

Then dimensionally reducing on the circle changes these terms to 

Sen 3 J d 5 x^7rM^a 1/2 (p) (l + |) R 5 + f ^ rfxy/ginMlro (l + | ) R 4 ■ (106) 

To go to the effective 'Einstein' frame, we absorb the radion perturbation by conformally 
transforming the metric to g\ B = Q 2 gAB- Although the graviton kinetic terms live in 
different dimensions, we can render the conformal rescalings of the two terms the same 
thanks to the fact that we are working entirely in the Gaussian-normal brane-fixed gauge. 
Since our starting metric satisfies g pp = l,g pp = 0, we can change the radial coordinate 
to p — j dp Q to ensure that the bulk term changes under the conformal transformation 

as dp ^fgl i? 5 = dp' ^Jgl; Rf/Vl 2 , i.e. the same as the brane term. We bear in mind that 
we will need to change rc M coordinates in order to keep the reduced metric in Gaussian- 
normal brane fixed gauge. So to absorb the factors oc $ from in front of the graviton 
kinetic terms (106) we need fl 2 = 1 + * . This defines our radial gauge transformation by 

P = -fro dp(l + f ) + r o = P + \!r dp<& where we have chosen the integration constant to 
ensure that the brane remains fixed, at p = r . This transformation is of the form of general 
diffeomorphisms (61) with x p — \ ff dp Then to keep the perturbation h pp from appearing 
and ensure that we are still in the Gaussian-normal gauge, we need to pick x' p — ~d p x p i or 
therefore x^ — ~ \ ff dp Jf dpd^Q (where we can ignore the overbars in these expressions to 
the leading order). 
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But now, from these expressions we immediately see that the gauge transformations 
vanish on the brane! Hence, the gauge fixing on the brane is unaffected by the dimensional 
reduction, and when calculating the fields along the brane we can just use the expressions in 
the original coordinates, ignoring altogether their behavior in the bulk 2 . Thus we can use the 
conformal transformation = (1+ *)<7 M „ on the brane. Now after Fourier transforming the 
decomposition theorem formula (74), and replacing h 4 by X using Eq. (78), the full metric 
perturbation along the brane in the original frame is 

h^ = ^ u -k»k^ + X5^, (107) 

where on the brane \I/ = AX — 2ro£. After our conformal transformation to linear order in 
perturbation, we find that h^ u = + ^ r/^, or therefore, using the expressions for $, \l> on 
the brane to express them all in terms of the variables X and £, 

^^-(-^+GKK 

where of course £ is given by (99). This is the effective 'Einstein' frame field of matter 
rings on the wrapped brane, which may yet contain extra scalar-like contributions from the 
helicity-0 modes. To see how it behaves we now need to find out the explicit forms of Green's 
functions Qyi and Q x . 



5.4 Fields of Static Sources Along the Brane 

Since we are looking for the static solutions, obeying = (0, k), which implies that all the 
momentum space Fourier transforms are formally Q{k) = 2ir5(k ) Q(k), we can factorize out 
the energy Green's function, set k° = in the formulas of the previous section, and drop the 
integration / ^ from the Fourier integral. We then solve the Green's function equations (91) 
and (103) using - again, as in the construction of the shock waves - the standard technique 
of sewing together the solutions of the homogeneous equation on either side of the brane, 
regular in the center and far away. 

Off the brane, both equations (91) and (103) look the same. The solutions which are 
regular in the center and at infinity are 



g = { AI «\ k ;) < " £r » ; (10 9) 

— * 

where now k = \k\. We sew them together at the brane by using the boundary condition 
found with the Gaussian pillbox integration of Eq. (91), that now gives different results for 
the two Green's functions, 

A £ TT = + J_ Ag = ^l [k 2 + ^Mt _ (3(k 2 )]g x + — A , (no) 

Ml TT M 6 4 ' x M| L 4M 5 3 r K n M 6 4 ' V ; 



2 Forearmed with the knowledge that the fields are localized near the brane, which will be ensured by the 
asymptotic boundary conditions on Green's functions. 
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with (3(k 2 ) given in (101), and then setting the homogeneous solution to zero. This yields 
1 I (k P< )K (k P> ) 



6 

and 

^{kp^K^kpy) 



(111) 



x 



Ml k h {kr )K x {k^) + h(kr )K (k^- b ) + *jL[W + g< - /?(fc 2 )]/o(fcro)*o(*£5) ' 

d 5 u 

(112) 

where, just as in Eq. (47), 

TO, \ jlo(kp)K (k^), p<r ; 

I {k P< )K {k P> ) = \ ) W/w 1 6 nM \ . ( 113 ) 

In fact the Green's function Qyi is identical to the Green's function that appears in the shock 
wave solution (46). This should hardly come as a surprise, because it is obvious from Eq. 
(90) that it has the same form as the shock wave equation when the stress energy source is 
traceless. The only difference between these two cases is in the domain of Fourier integration, 
which for relativistic sources involves one integration fewer because of the Lorenz boost of 
the longitudinal direction. 

Since we wish to explore the fields along the brane, we set p = r in the Green's functions 
(111) and (112). At distances £ ^> r , where we can limit our attention to axisymmetric 
configurations, for which the solutions (111) and (112) are valid, kr^ ~ r^/i <C 1, we can 
always replace Iq and I\ by their small argument expansion, and approximate them by 
Iq — > 1, I± — > 0. We must be more cautious with K and K\ since they depend on the deficit 
angle, and for near critical branes kr / (1 — b) may be big even when kr <C 1. Thus we will 
need to consider the limits b < 1 and b — > 1 separately. We have already encountered this 
previously, in the example with shock waves. With this in mind, we see that the consistent 
approximation of the Green's functions (111) and (112) is 

1 1 

1 1 

The TT-tensor sector is particularly simple to understand from (112). As with shock 

i ill- c i i • ill -7 M £ K l( k T^h) 

waves, we see that the behavior of the theory is governed by the ratio k c = ,, 3 „ ,. r Q \ , 
which controls the denominator of Qtt at low momenta. For k < k c , or distances £ > l/k c , 
the dominant contributions always come from oc „ ,,. 1 >~a b \ ■ In this limit the long range fields 

"■0(K 1 _ b ) 

manifestly reveal the extra dimensions, since the scaling of the potentials will not be AD. In 
the generic sub-critical cases, b < 1 and for momenta < l/r we can replace the Bessel 
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functions K v by their small argument expansion along the brane, which yields the crossover 
scale identical to Eq. (52), r 2 c (k) = 27r (i-b)M i ^ n 2 fcr^ ' beyond this distance, the theory 
behaves as a 6D gravity on a cone. 

In the near-critical limit, however, b — > 1, at intermediate momenta we should instead 
approximate K v by their large argument expansion, as in the shock wave analysis of the near- 
critical models building up to Eq. (56). In this case the crossover scale is r c = jfi, where 

6 

gravity first changes into a 5D theory, because for near-critical tensions one bulk dimension 
is efficiently compactified on a circle, as discussed recently in [59], and in our shock wave 
analysis. Eventually, the circle opens up and the theory again turns into a 6D gravity on a 
cone. As a result we see that the crossover scale, beyond which the theory will not look 4D, 
is exactly the same as in Eq. (59), which was already revealed by the shock waves. We again 
confirm that the see-saw mechanism of [29] is realized only in the near-critical limit, but the 
crucial dynamics which manufactures the see-saw scale is the compactification of one bulk 
dimension induced by a near-critical brane [59]. 

What remains is to check precisely what kind of AD theory we have below the crossover 
scale r c , in the regime of distances tq <C £ <C r c . For TT-tensors, by previous discussion 

we can neglect the terms in the denominator ~ ,, Vt^, . The Green's function reduces to 

Gtt — > M l k 2 ■ The non-relativistic matter sources on the wrapped brane are described by 
stress energy of Eq. (35), which Fourier-transforms to r a 5 = — /i diag(l, 0, 0, 0, 1). So the 
leading order solution for the TT-tensor field is 



2/x 



ji — v — , 



Next we need to evaluate the Fourier integrals 7 M „(:r) = / -^^l^u{k)e lk ' x . Since a Fourier 
integral picks up the dominant contribution from the momenta k ~ in the AD regime, 

when \x\ < r c , we can approximate 7 M „(/c) by (116). The remaining integral is easy for 7%, 
and yields 3 7% = 3n ^^ in this regime. Recalling that fi is the mass per unit length of string 

(see Eq. (31)), so that ji = and using as before Gauss law to trade the 5D Planck scale 



M 5 for the AD one, Mf = we finally obtain 

M 

3TrMi\x\ 

To find 7 J, fc, we don't need to do the Fourier integrals directly. Instead, using Lorentz 

invariance, we see from (116) that 7^ = j^(^ j k + £2^=5^- Then since 7^ is a TT-tensor, 

and 7 J, = 0, using transversality djrfk = 0, we get C 2 = 1. Further vanishing trace 7 M M = 
implies ^ k k = — 7°o, which by comparing with Eq. (117) sets C\ = — X2 ^ M 'x ■ Thus we find 

^ = -™(^ + ^)- (118) 



7°o(^) = T7rW- ( 117 ) 
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It is convenient to rewrite the solutions by introducing the effective AD Newton's constant 
Gn eff = %^ M 2 ■ Then, in the leading order for |x| < r c , 



7°o(x) 



3° Neff W 



-- G 



M 

1" ' 

M 



3~ Neff W\ 



5 j k + 



x J x k 

X 2 



(119) 



If we were to ignore the scalar field X for the moment, the Newtonian potential Vn = — h°o/2 
would be Vn = — § G^eff which is a factor of 4/3 larger than the usual formula of General 
Relativity. This is a manifestation of the Iwasaki-van Dam-Veltman-Zakharov discontinuity 
[17], which signals the presence of the helicity-0 modes in the graviton multiplet. The factor 
4/3 enhancement is precisely what one expects based on other examples [15], and in this 
case shows that the extra helicity-0 modes in the spin-2 multiplet mediate attractive force, 
just as the helicity-2 modes. Thus, they are not ghosts. 

The scalar Green's function Q x is considerably more involved because of the function 
f3(k 2 ) in the denominator. It behaves differently for generic sub-critical and for near-critical 
branes, and so we need to deal with it with some care. First off, we note that the terms 



k 2 + jjT^t — I3{k 2 ) factorize as 



4M|r 



k 2 + 
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4M, 3 



5^0 



-(3{k 2 ) 



r 2 k 2 



l-b + rlk 2 



1 - 



3(1 -b)M; 



46M 6 4 r 



k 2 + 



3bM* 



V4M, 3 



5?"0 



1-6 
~2r 2 T 



(120) 



After substituting this in Eq. (115) we find that Q x reduces to 



1 



1 



k 



Ml 
7T 



rik 2 



Ml 1-b+r^k 



TT71 



4bM*r 



k 2 + 



4M?r 



1-1 
2^ 



121) 



Now, we will generically require r <C M|/M| in order to have the crossover scale for gravity 
be much larger than the radius of the compact dimension, r c ^> r . Otherwise, we can never 
get the theory to behave as a AD gravity. 

For generic sub-critical branes, b < 1, at scales k <C l/r the scalar Green's function 
(121) is approximated by 



g 1 rolnPffi 



;i22) 



where we have approximated .,V \ with the small argument expansion. Hence at very 
large distances, or low momenta, the scalar Green's function depends on the momentum 

, 1-2(1-6), 

only logarithmically, Q x — , implying that the configuration space solution for the 

scalar field X depends on the distance as X ~ M/\x\ 3 - i.e. as a field in QD. However, 
because the momentum-dependent terms in the denominator of (/ x are negative, as we move 
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inwards towards the source, the scalar force grows, and the scalar eventually becomes strongly 
coupled as k approaches the pole of (122). This happens when k ~ fc*, where 



2(1 - b)M ( 



(123) 



But this is essentially the same as the crossover momentum in the TT-tensor sector, as we 
can see by comparing k* with k c given in Eq. (50), after recalling that Mj ~ M|r . At 
first sight, it looks a little surprising that for resolved codimension-2 branes the strongly 
coupled scalars appear even around the vacuum, in contrast to the codimension-1 brane 
induced gravity models, as discussed in [19, 20, 21, 22, 23]. However the AD vacua aren't 
really vacuum solutions from the point of view of the full brane worldvolume theory because 
they contain the flux of the axion E and hence source the brane bending term £ even in the 
absence of localized matter sources, as seen in Eq. (99). This in turn triggers the onset of 
strong coupling. 

Thus we see that for sub-critical branes with b < 1, inside the regime of length scales 
£ < r* where 

9 Mj , ,2(1-6), , x 

the TT-tensor is AD, but with the wrong tensor structure, as we have discussed following Eq. 
(119), however the scalar sector is in fact strongly coupled. This makes the perturbation 
theory around the vacuum on sub-critical branes completely unreliable inside the regime 
where gravity might be AD. The negative signs of the momentum-dependent terms in the 
denominator of (122) might naively suggest presence of ghosts in this regime, but at this level 
of the approximation we cannot conclude that decisively. One exception to the pathological 
behavior of the theory is the limit of relativistic sources, for which = = 0, so 
that the scalars are never sourced and the tensor sector reduces to the AD Aichelburg-Sexl 
solution below the crossover scale. One needs some means for exploring the scalar sector 
beyond linearized theory, perhaps along the lines of [68], or by seeking exact solutions, 
before passing on the final verdict on the sub-critical brane models. At the level of linearized 
perturbation theory the gravitational effects below the crossover scale are not calculable. 

In the near-critical limit the scalar sector behaves very differently. At distances I ^> tq, 
where kr ~ ^ C 1 the argument of the Bessel functions K v will still be very large when 
b — > 1. We have already noted this in the discussion of shock waves. Hence for distances 
below the crossover scale we can use the large argument expansion for K u 's, which yields 
K\ — > Kq, as long as r$k ^> 1 — b. Moreover, we can immediately neglect the negative terms 
in the denominator of Q x in Eq. (121) because ^^I^ 5 = (1 — b)^ < 1. The scalar Green's 
function is therefore ^ 

^ x = T7I , ; Wl rgfea — • ( 125 ) 



Mj l-b+rfik 2 V ~ 4r c ro' 

where r c = M|/Mg is the TT-tensor crossover scale, as per Eq. (56). Now, it is convenient 
to define the scale r vac by 



„2 



b=-?-, (126) 



'"vac 
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In the limit b — > 1 we have r vac = Vl — *C and so the scale r vac is always smaller 
than the size of the conical throat surrounding the brane. This is why the scale r vac may 
compete with r c for the control over the scalar sector. Both r vac and r c are smaller than the 
size of the throat, and the details of sub-crossover dynamics depend on their ratio. First, 
note that for momenta k such that r vac /c > 1, we can expand 1 _ b ° hr 2 fc2 — 1 + C(^zp) and 
keep only the first term. The scalar Green's function in this limit at momenta r k < 1 — b 
reduces to 

ex = WiT^r- (127) 

Thus for momenta k > max(l/r c , l/r vac ), the theory looks 4D, but the scalar has a mass 
term, m\ = ^r^- Therefore its long range effects have a Yukawa suppression oc exp(— m x r). 
Hence in the range of scales ^ < £ < min(r c ,r vac ) the scalar field exchange will lead to 
negligible, exponentially suppressed effects compared to the TT-gravitons. This mass term 
is completely analogous to a brane-localized mass term of a bulk scalar. In such models the 
scalar is repelled from the brane, since it prefers to reside in the region of space where its 
inertia is minimized. This resembles Meissner effect in superconductivity, where the mass 
term pushes the magnetic field outside of the superconducting medium. In our case however, 
the scalar is also pulled back to the brane by the brane-localized kinetic terms, and the result 
is the Yukawa suppression alone, which is manifest in (127) in the same coupling 1/Mf as 
for the TT-tensors below the crossover scale, e.g. in Eq. (116). Similar phenomena were also 
studied recently for gauge fields in gaugephobic models [69]. 

Although X may be Yukawa-suppressed, so its direct long range forces will be small 
compared to those arising from the TT-tensor sector, the brane bending term £ may yet 
yield a long range tail, as is clear from Eq. (99). We see that £ is proportional to k 2 X, 
and this will in general lead to long range effects that may dominate the Yukawa-suppressed 
effect from direct X-exchange. It turns out that the ^-induced effects will be parametrically 
smaller than the dominant TT-tensor fields, as we will show below. 

At distances £ > r c , the linear term in k in the denominator of (127) dominates the 
momentum transfer in scalar exchange, modifying the theory to a 5D one, and eventually, 
after resuming the corrections which we neglected in setting K\/Kq ~ 1, one would see how 
the theory changes to 6D. When r vac > r c the story of the scalar dynamics ends here. 

However, when r vac < r c , in the regime of scales r vac < I < r c below the crossover scale, 

r 2 k 2 r 2 k 2 

we can approximate the fraction in the denominator of (125) by 1 _ b + r 2 k 2 — ~f^- Thus, 
Qx 1 ~ k + jYZ^M^k 2 {k 2 + 47^-)- Then, since the mass term m\ is larger than the momenta 

in this regime, because k 2 < 1/V^ ac so that ^ < j^P 1 = (1 — 6)^ < 1, the Green's function 
is approximately 

Mq k + ^r^k 2 

What happened here is that due to the momentum dependence of the effective coupling, the 
scalar field mass term got overcompensated and the field became essentially massless all the 
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way out to a new scalar crossover scale, given by (using Eq. (126)) 



^ = z— = zrV ^ 

4 r 4 1 — o 

which is essentially equal to the length of the throat inside which bulk gravity looks 5D. 
Below TZ C , the scalar propagator reduces to Q x 3M 4^° k 2 < M l k 2 , because r c < ^ in the 
near-critical limit. Thus, when r vac > r c , the scalar remains massless all the way out to 1Z C , 
but below the crossover scale it couples more weakly than the TT-tensor. 

We can now put together the various contributions to the fields of a static massive ring on 
the brane at distances below the crossover scale. The scalar Green's function for near-critical 
branes below the crossover scale is well approximated by 



w .i /,2_,„2 • r c < r vac, l c < k , 

3M|r F ' rvac < r c j; < k < 



The scalar field X sourced by a ring of mass with Fourier-transformed stress energy tensor 
T a b = — /xdiag(l, 0, 0, 0, 1) is then, using the solution (104) for this range of scales, and again 
replacing /x = and Mf = as in the earlier discussion of TT-tensors, 




*^ '"vacj „ *^ k , 



A - 1 2A4 (l-6)r c 1 „ < J_ < fc < J_ 

^vac 



while the brane bending, normalized to ro, is, from (99), 

t [ r ° r " M (lO 1 \ < ± < f c 

? J ^v 2 ac 6Mj { k 2 k 2 +m 2 h ' c ^ ' vac ' r c ^ ^ > 



(132) 

r (i^K i + constant r < r c , — < A; < — . 

The constant in the brane bending transforms to a gauge- dependent contact term oc 5^ (x) , 
which is ultralocal and hence we can neglect it altogether. The configuration space solutions 4 
for the scalars X and £ are therefore 

M p -mx\x\ \r\ < r < r 

!2nM 2 \x\ e ' Fl <~ r c<~ 'vac , 

(l-fe)^c M _ /|^| <r 

and 



'vac -i«4 



vac 5 



(134) 

L r 12^Mjp| ' 'vac ^ K| ^ 'c • 

Now, in the case when r c < r vac , the scalar field X will have a very short range due 
to its mass. Indeed, since m~^ c ~ r r c is the geometric mean between the brane radius 
and the crossover scale, m~^ c will be quite small. In fact, if the brane size is set by the 



4 Which wc find using / ( ffi„ = ^ 
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AD Planck length and the crossover scale by the current horizon size, r ~ 1/M 4 , and 
r c ~ 1/H respectively, m^ 1 will be near the table top bounds of about 0.1mm. Thus X is 
de facto decoupled at the scales where we normally probe 4D gravity. On the other hand, 
the brane bending can still simulate a long range effect which scales as l/\x\ due to the 
momentum dependence of the couplings. However, the strength of its contributions is a 
factor ^ = ^ i down compared to the TT-tensor contributions in Eq. (119). In 

''"vac ?*vac ?*vac v y 

the case when r vac < r c , the scalars remain essentially massless and source long range fields. 
However, their strength is still suppressed, this time by ^~^ rc <C 1, by definition of the 
near-critical limit. While the scalar effects become more important as this factor becomes 
smaller, so do the strong coupling effects which we cannot neglect away from the near-critical 
limit. Indeed, as ( - 1 ~^ rc — > 1, where the brane bending effects could compete with the long 
range fields of TT tensors, the strong coupling effects become very important, and we lose 
calculability in linearized perturbation theory below the crossover scale. 

Hence in either case the long range effects of the scalars in linearized perturbation theory 
around the near-critical vacuum cannot compete with the TT-tensors, as long as the pertur- 
bative treatment is valid. We can therefore ignore the scalars in the formula (108). The 
long range fields are determined by 7^ alone, which are given in Eq. (119). Clearly, these 
solutions do not look like 4D General Relativity. Instead, in the leading order they really 
mimic a scalar tensor theory, due to the presence of the helicity-0 contributions. To see this, 
we can compare the solutions (119) with the spherically symmetric linearized solution in the 
PPN approximation [64]. For the purpose of this comparison, we define G^ = fGVe//) so 
that 
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= 2Gnt-^t 


\x\ 






l 3 k 


- ~G N - 



\4 1 ~ MxH k 

J\ 5k ~2 N ^xf~' ( ) 

Comparing to [64] (taking into account overall sign difference reflecting our conventions) we 
find that the solution (135) mimics a Brans-Dicke theory with 7 = |, or therefore, with the 
Brans-Dicke uo parameter equal to zero, where the parameter u is defined in the usual way 
in the Brans-Dicke action as Sbd = / d A x^fg(QR — uj(VQ) 2 This theory is in conflict 
with observational data as it stands, and hence is not a realistic description of our Universe. 

However, it is possible that non-linear effects, which we have neglected throughout this 
work, could play a role here. Indeed, even in conventional General Relativity, nonlinearities 
start to show up at distances comparable with the gravitational radii of the sources. In brane 
induced gravity, and other frameworks that strive to modify gravity in the IR, the nonlinear 
corrections will become important even sooner than the nonlinearities in General Relativity 
[18, 19, 20, 21, 22, 23]. Our results are found in linearized perturbation theory around a 
near-critical vacuum, and so it is conceivable (although not certain) that different strong 
coupling effects at higher orders in perturbation theory might improve the behavior of long 
range fields of masses far from the source but well below the crossover scale. Alternatively, 
it is possible that additional tweaks of the bulk theory, for example by curving the bulk 
locally, could change how the helicity-0 mode couples to brane matter. Such methods are 
widely used in the construction of string landscapes [3], and might also be useful here. 
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We stress again that regardless of the ratio of r vac /r c the scalar couplings remain consis- 
tently weak for static sources on near-critical branes with b — » 1. The static solutions are 
governed by Euclidean momenta, and so they are always finite in linearized theory around 
near-critical vacua, in contrast to the fields of masses which perturb generic sub-critical 
vacua. Thus linearized perturbation theory remains under control on near-critical vacua. 
This as we already mentioned does not guarantee that the linearized solutions will remain 
dominant close to the masses that source the fields. Further study of these issues, to ascertain 
how likely it is that nonlinear phenomena may yield any additional screening of helicity-0 
modes, hence seems warranted. The near-critical vacua at the very least provide us with a 
controllable new arena where such phenomena could be studied. 

The near-critical propagating solutions, on the other hand, also remain under control, 
albeit their dynamics is more subtle. The issue is that there may be new poles at Lorentzian 
momenta k 2 = —k 2 = —^r, in the solutions for £ and X, given respectively by Eqs. (99), 

(104) and (112), and their approximations below the crossover scales. Clearly, these poles 
never play any role for static configurations, which are controlled by Euclidean momenta 
for which the fields remain finite. However one may worry if new infinities could plague the 
linearized theory once a weak time dependence is allowed. Now, on a generic sub-critical 
brane the poles at k 2 = — &^ ac reside practically on the 4D cutoff ~ 1/Vo, deep in the strongly 
coupled regime of the sub-critical linearized theory and hence can be completely ignored at 
large distances as a UV mirage. On the other hand, in the near-critical limit b — > 1 the 
scale /c vac becomes very low, and may be much lower than the crossover scale l/r c . However, 
as we will now explain, in this case these poles are harmless because they merely point to 
a breakdown of a gauge fixing of the linearized theory, which is fixed by an interchange 
of scalar modes. The easiest way to see this is to reconsider the formula for the £ field, 
Eq. (99), and the boundary condition for AX', Eq. (100), when Lorentzian momentum 
is k 2 = —k 2 = —^r. Then, Eq. (99) indicates that £ may diverge; however since the 
brane bending must be bounded, this equation must be reinterpreted by multiplying it by 
1 — b + r\k 2 before setting k 2 = —k 2 ac , and then demanding that £ is finite, which means 
that the RHS must vanish identically. This fixes X on the brane to 



X 



2r, 



(136) 

1.2— _U2 



while now £ is not determined by this boundary condition. Since X is now fixed by the 
matter source, it cannot be freely chosen to satisfy the boundary condition for AX' in (100). 
This is what fixes £: indeed, using (99) and (101) we can rewrite (100) as 



AX' 



fc2=-fe T 2 ac M| 



bM*r \ £ 1-b 



Ml ) 2rg \-b- &J g;° 6M| 6Mf 
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so that for given sources and for X fixed by (136) when k 2 = — ^ we can choose £ to satisfy 

(137). Therefore the Lorentzian poles at k 2 = —k 2 ac are spurionic, and can be ignored in 
the linearized perturbation theory on near-critical branes. The linearized theory around the 
vacuum on such generic near-critical backgrounds remains under control, without dramatic 
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instabilities, and mimicking Brans-Dicke theory with u — 0, although we stress again that 
one has to reevaluate it against the higher order corrections closer to the source, because of 
the issues related to the Vainshtein scale and strong coupling at higher orders. 

6 Conclusions 

Our main result is the observation that properly regulated brane induced gravity theories, 
which yield calculable long distance gravitational fields, really behave as a semiclassical 
landscape of vacua. Treating gravity classically, these models can impersonate AD worlds 
between the UV cutoff that resolves the core of the gravity-localizing defect and the crossover 
scale, if their Scherk-Schwarz sector is carefully tuned to yield AD Minkowski vacua. While 
we have explicitly worked with codimension-2 defects, we feel that similar conclusions should 
extend to any setups with codimension > 2. In the case of codimension-2, we see that the 
brane vacua can remain completely flat, while the bulk readjusts to absorb the tension as 
deficit angle. Separating them there will be mismatched configurations, where the tensional 
pressure in the compact direction may not be precisely cancelled, which are either nonsta- 
tionary, or never approximate AD behavior. Static supercritical branes must be singular, just 
like supercritical local strings in AD [54]. However in this case there should exist nonsingular 
inflating solutions, that should really be the correct vacua for supercritical branes. 

We also find that the cosmological constant problem changes its guise rather dramatically 
in codimension-2 setups. Although the theory has many AD flat vacua with Minkowski 
metric, the effective AD Planck scale, and for sub-critical tensions, the crossover scale out to 
which gravity looks AD, are very sensitive to the brane tension. These parameters depend 
on the tension directly, as does the sub-critical crossover scale (see Eq. (59)) or through 
the compactification radius of the wrapped 4-brane, as does the effective AD Planck scale. 
So if the theory is to mimic a weak AD gravity over a large spatial region, which looks flat 
and static, for a given brane and bulk Planck scales M 5 and M 6 and a fixed brane radius 
r one must tune the brane tension precisely. Changes of brane tension will generically 
lead to proportional changes of the crossover scale and AD gravitational coupling, requiring 
appropriate retunings. However in the near-critical limit, where the bulk compactifies to a 
cone so that AD gravity first changes to a 5D one, the crossover scale saturates at a value 
completely independent of the tension, r c = M|/M|. To get this to be of the order of the 
current horizon size, Hq 1 ~ 10 28 cm, and ensure that the rings of matter on the brane look 
pointlike at energies below a TeV, one needs M 6 >, TeV, and therefore M 5 <, 10 19 GeV. 
Further to reproduce the AD Newton's constant below the crossover scale, one needs to pick 
r £j 10~ 19 GeV _1 . These numbers are rather curious. Clearly, one must tune the theory to 
make sure that the brane is very thin, with many hidden sector fields, to get such hierarchies. 
We don't have much to add to the discussion of how to do it, but merely note that this is in 
line with the current philosophy of the brane induced gravity models [15]. While some work 
along the lines of embedding the theory into string theory, that would provide the framework 
for its full UV completion, has been pursued [60, 61, 62], the existing constructions are really 
semiclassical models where gravity is treated classically. Adopting this possibility, one may 
at least explore the low energy consequences of such models. 
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Once we consider gravity of localized static masses on the brane, we find that for sub- 
critical tensions there is a strong coupling scale in the vacuum itself, which is essentially the 
same as the crossover scale. Thus in the regime where gravity may appear 4D, linearized 
perturbation theory does not apply because the scalar modes are strongly coupled. However 
in the near-critical limit the linearized perturbation theory around the vacuum is under 
control. We find that it does not contain instabilities, and generates long range fields which 
are, to the leading order, the same as in the Brans-Dicke theory with uo = 0. This would be 
in conflict with tests of gravity in Solar System and beyond. However, throughout this work 
we have taken the bulk to be locally flat. Adding bulk fields and curvature will increase the 
diversity of possible solutions and might yield additional effects screening the helicity-0 mode 
in the spin-2 sector, perhaps similarly to what happens in string landscape constructions [3] . 
Moreover, the nonlinear corrections might play some role too as advocated in [18, 19]. 

This discussion shows that many questions remain open. Clearly, the most interesting 
phenomenological issues concern the stabilization or decoupling of the helicity-0 modes. It 
is interesting to explore the theory for dangerous instabilities beyond the linear order. It 
also remains to see if the brane induced gravity models can be consistently derived from 
microscopic models that admit plausible UV completions, generating M4 ^> Mq. In other 
words, do such landscapes even exist beyond a classical action, that one can write for them? 
Finally, it would be interesting to see if the new guise of the cosmological constant problem, 
and in particular the serendipitous insensitivity of the crossover scale from brane tension 
in the near-critical limit, might yield some new avenues for going around the venerated 
Weinberg no-go theorem [1], that is still the main obstruction to having a mechanism for 
protecting a small 4D vacuum curvature from quantum radiative corrections. These issues 
are outside of the scope of the present work, and we can only hope that our sketch of such 
a diverse new landscape with a different manifestation of the vacuum energy problem may 
motivate the search for their resolutions. 
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